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Abstract. We review recent simulation studies of interfaces between immis- 
cible liomopolymer phases. Special emphasis is given to the presentation of 
efficient simulation techniques and powerful methods of data analysis, such 
as the analysis of capillary wave spectra. Possible reasons for polymer in- 
compatibility and ways to relate model dependent interaction parameters to 
an effective Flory Huggins parameter x a- re discussed. Various interfaces are 
then considered and characterised with respect to their microscopic structure 
and thermodynamic properties. In particular, interfaces between homopoly- 
mers of equal or disparate stiffness are studied, interfaces containing diblock 
copolymers, and interfaces confined in thin films. The results are related 
to the phase behaviour of ternary homopolymer/copolymer systems, and to 
wetting transitions in thin films. 



1 Introduction 

Blending chemically different polymers is a cheap and relatively straightfor- 
ward way of creating new materials, and polymeric alloys are therefore indus- 
trially and technologically omnipresent. Prominent examples are, e.g., rubber 
toughened plastics. Their widespread use notwithstanding, polymer mixtures 
are seldom homogeneous at temperatures of practical interest. Any slight in- 
compatibility of the monomers, as is usually present between different organic 
molecules, is amplified by the number of monomers in the macromolecule, and 
cannot be balanced by the entropy of mixing for typical molecular weights 
[0 . On a mesoscopic scale, such materials consist of numerous microdroplets 
of one phase, which are finely dispersed in the other phase. The material 
properties thus depend sensitively on the structure and the properties of the 
interfaces between different phases. Correlations between interfacial and bulk 
properties are present on various length scales. On the one hand, the local 
interfacial structure - the conformation of polymers, the interfacial width, 
which is closely related to the number of entanglements between polymers of 
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different type, enrichment of chain ends or solvent at the interface etc has a 
fundamental influence on the mechanical stability of the alloy. On the other 
hand, the morphology of the blend at given conditions of preparation (e.g., 
given stirring rate during mechanical mixing) is basically determined by the 
interfacial tension: An old theoretical argument due to Taylor]^] balances the 
viscous stress (caused by stirring) and the interfacial tension, and predicts 
that the droplet size should be directly proportional to the latter. This law 
is indeed found experimentally 

A huge number of interesting questions are connected with the general sub- 
ject of polymer interfaces. For example, interfacial properties can be tailored 
by adding a small amount of a third substance to the blendQ. In particular, 
copolymeric surfactants containing both types of monomers are often used 
as effective compatibiliser. Their effect on the morphology of the blend is 
twofold. First, they reduce directly the interfacial tension: Being compati- 
ble with both components, they aggregate at interfaces, thereby reducing the 
number of direct contacts between homopolymers of different typej^]. Second, 
recent experimental^ and theoretical[Q studies indicate that their presence 
at the droplet surface prevents the coalescence of the droplets brought into col- 
lision in the course of mixing. The effect can be related to the two-dimensional 
compression modulus of the copolymer filmQ]. Since droplets can break up, 
but do not merge any more, one obtains a particularly fine dispersion. In 
addition, copolymers improve the mechanical properties of interfaces. The 
interfacial width increases, and likewise the number of entanglements, they 
adhesive attraction and the fracture toughness At high enough copoly- 
mer concentrations, additional copolymer rich phases emerge which display a 
variety of structures ordered on a mesoscopic scale || |nj, [□], [T^]. These meso- 
scopically structured materials promise to possess unique and useful materials 
properties 0, 

Another complex of important problems refers to confined interfaces, i.e., 
interfaces interacting with one or two surfaces. The presence of a surface 
influences the interface on all length scales - the local structure is affected 
as well as long wavelength fluctuations of the interface position. Depending 
on the interactions of the surface with the different components of the blend, 
the interface may bind or unbind, giving rise to a whole diversity of wetting 
phenomena [u| 0. 

From the point of view of basic science, inhomogeneous polymer systems 
are interesting because of the different length scales involved (polymer gyra- 
tion radius vs. monomer size), and because of the additional conformational 
degrees of freedom of polymers as opposed to smaller molecules. Interest- 
ingly, these apparent complications have partly the effect of simplifying the 
physics: Since polymers interact with so many other polymers, microscopic 
details of the chemical structure of monomers wash out to a large extent, and 
can be absorbed quite successfully into a few number of effective parameters. 
Furthermore, the effective interaction range - roughly the polymer gyration 
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radius - is very large, and the region in which critical concentration fluc- 
tuations become important is extremely small as a consequence (Ginzburg 
criterion, see Rcf. fTifl). Polymer blends are thus unusually well described by 
mean field type theories. 

On the other hand, the treatment on the mean field level in itself is al- 
ready very involved, especially if one attempts to account for local correla- 
tions, and is interested in local structures. A number of mean field type 
approaches have been established in the past years, which differ by level of 
coarse-graining and by the type of questions they address. Among these, we 
quote Landau-Ginzburg [[l9[ |2(| and scaling approaches Q, which coarse- grain 
over the microscopic and to some extent even over the chain conformational 
structure; Self-consistent field theories ^2|, |2^, ^ and density functional 
theories p5l , which treat chains as random walks in a mean field environ- 
ment, mostly ignoring the monomer structure and local chain correlations; 
lattice-based theories such as the famous Flory-Huggins theory |2(| |2^] and 
subsequent more refined extensions (e.g, p§, |2^, 30 ); and finally Schweizer's 
and Curro's P-RISM theory which incorporates the local liquid structure 
into a theory of polymer melts, using concepts from integral equation theories 
for simple liquids. The highly coarse-grained theories have the advantage of 
relative simplicity, thus allowing in many cases for an analytical treatment. 
However, the information they can provide on local structure properties is, 
a priori, very limited. Furthermore, they require a number of "effective" 
parameters as input, which cannot be determined from microscopic param- 
eters within the theory. On the other hand, the more microscopic details 
are incorporated into a mean field theory, the more involved the treatment 
gets, and the more additional approximations have to be made in order to 
make it tractable at all. At the lowest level of coarse-graining, for example, 
the P-RISM equations totally neglect chain end effects, but their truly self- 
consistent solution in a one component melt nonetheless requires a series of 
single chain Monte Carlo simulations |H[ Ej2]. 

The universal aspect of the physics of polymeric alloys, and the close re- 
lationship between the local structure of interfaces and the global material 
properties, make them particularly suited for computer simulations [B^ 
These provide simultaneously a detailed microscopic picture of the interfacial 
structure, and information on the thermodynamics of the interfaces. When 
compared to experiments, they serve as a test of the microscopic model which 
has been used. When compared to theories, they serve as a test of the theory, 
within a well-defined microscopic model. In addition, they may provide struc- 
tural information which may not yet be accessible experimentally or theoreti- 
cally. They can thus contribute substantially towards a deeper understanding 
of the connections between the microscopic parameters, the microscopic struc- 
ture and the macroscopic properties of a material. 

Obviously, there are limitations. With the present computational re- 
sources, full calculations of polymer melts in atomistic detail are currently 
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far beyond reach. It is thus necessary to take again advantage of the univer- 
sality ideaE |2?j, [35|, and represent real polymeric system by idealised polymer 
models. In this spirit, a number of chemical monomer units are mapped onto 
one effective monomer with a much simpler structure. Even within such a 
simplified model, only melts of polymers of rather modest chain lengths can 
currently be equilibrated and subsequently studied. Depending on the ques- 
tions that one wishes to study, models of different levels of idealisation have 
to be chosen. The careful choice of a suitable model is thus crucial for the 
success of an investigation. We will discuss this important point in more detail 
in Sec. 3. The reader interested in a general overview over the use of computer 
simulations in polymer science is referred to, e.g., the set of excellent reviews 
in Ref.|§. 

In this contribution, we will review some recent simulation studies of poly- 
mer interfaces in polymer blends. We shall restrict ourselves to interfaces at 
thermodynamical equilibrium, and to studies of static properties. Our survey 
will focus on the insight that simulation studies can provide into the local 
structure of interfaces, and the implications for the global thermodynamics of 
the systems. We will start with the closely related issue of the general rela- 
tionship between local correlations and thermodynamic miscibility in binary 
blends. Then, we will discuss interphase boundaries in binary homopolymer 
blends. The effect of adding diblock copolymers on the interface and on the 
phase diagram is considered. Finally the behaviour in a thin film is exam- 
ined, with special emphasis on the effect of the interfacial fluctuations on the 
measured profiles. 

2 Polymer incompatibility and Flory Huggins 
parameter 

We begin with a discussion of polymer incompatibility: Immiscibility in poly- 
mer blends can be caused by several factors. First, one has usually a direct 
relative repulsion between monomers of different type. In nonpolar molecules, 
for example, the van der Waals attraction between monomers i and j is pro- 
portional to the product of their polarisabilities otiOtj. Thus the interaction 
between unlike monomers is smaller than the arithmetic mean of the in- 
teractions between like monomers - which gives rise to a relative repulsion 
cx {&A — <5b) 2 , where 5, oc on is the Hildebrand solubility parameter |Q. This 
enthalpic incompatibility is inversely proportional to the temperature T. It 
may be supplemented by entropic effects: If the monomers have different 
shapes, like monomers tend to pack more efficiently than unlike monomers. 
This effect has been studied in detail within a lattice model by Freed and 
coworkers p8[ . Similarly, stiffness disparities tend to favor demixing, as has 
been shown within P-RISM theory j37j: When mixed with stiff polymers, the 
more flexible polymers loose conformational entropy, and even though the 
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stiff polymers win entropy in return, the net effect turns out to be negative. 
Fredrickson, Liu, and Bates pointed out that phase separation in blends with 
components of different flexibility or architecture is also promoted by long 
range composition correlations j38|. The entropic contribution to the incom- 
patibility of polymers has no direct temperature dependence. An indirect 
temperature dependence may enter through the chain stiffness. The sum of 
entropic and enthalpic contributions will thus generally lead to a complicated 
temperature behaviour |2£}|. 

The incompatibility of polymers in a binary A/B-blend is often described 
in terms of a single Flory-Huggins parameter \- It has originally been derived 
from a simple lattice model on the base of three assumptions H: 

(i) The distribution of polymer conformations does not depend on the com- 
position of the blend. 

(ii) Composition correlations are neglected. 

(iii) All monomers have equal size (one lattice site per monomer), and the 
melt is incompressible (i.e., the lattice is fully occupied by monomers). 

The approximation (ii) implies, first, that nonlocal correlations induced by 
the chain connectivity (the "correlation hole" ) are ignored, and second, that 
short range composition correlations related to local demixing tendencies are 
disregarded ("random mixing"). One obtains the free energy of mixing per 
site 

F FH /(k B T) = ^ ln(^) + ^f- ln(0 B ) + xfate, (1) 

where 4n is the volume fraction and iVj the chain length of component i, and 
the last term describes the loss of enthalpy upon mixing. Assuming that 
neighbour monomers i and j interact with the interaction energy ejj, the 
Flory-Huggins parameter \ is given by 

z-2 ( e AA + e BB \ ,„s 

X = J B Y[ €AB 2 J' (2) 

with the coordination number z of the lattice. Equation (|^) takes into account 
that the interaction of a monomer with its two neighbours along the same 
chain should not contribute to the energy of mixing. The free energy (Q) 
describes usual demixing behaviour, with a miscibility gap at x parameter 
values larger than 

In the one phase region, the \ parameter can be determined, e.g., from the 
small angle structure factorpg]. As long as \ is used as an adjustable, heuristic 
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parameter, this simple model is found to perform extremely well in comparison 
with experiments. The value of x itself is theoretically less accessible. For ex- 
ample, neutron scattering data reveal in apparent contradiction with Eqm(0) 
that x mav depend significantly on the composition of the blend [ll], |42| . 

However, this apparent failure is not surprising, since the crudest version of 
the Flory theory disregards a number of effects which influence real polymeric 
fluids. In particular, the relative repulsion between polymers of different type 
depends on the number and the character of contacts between polymers. One 
can distinguish between different factors. 

First, the total volume (at given constant pressure) or the pressure (at 
given constant volume) of the mixture depends on the composition of the 
blend. This effect is commonly referred to as the "equation-of-state" ef- 
fect. Experiments are usually conducted under constant pressure conditions, 
whereas simulations and analytic calculations often use a simulation box of 
fixed volume. Since the thermodynamics of mixtures is conceptually simpler 
at constant pressure p, we shall mostly discuss the NpT ensemble in this 
section. We will relate our conclusions to constant-volume simulations later. 
At a given pressure, the different components will generally have different 
monomer densities g* A and g* B in the pure states, and the volume per monomer 
in the mixture is to lowest approximation the appropriately weighted average 
of l/g* A and l/g* B . In addition, the incompatibility of unlike chains often 
leads to an "excess volume on mixing" v cxc . As we shall see, the composition 
dependence of v cxc is roughly parabolic, i> cxc ~ vp(l — p). All taken together, 
the density g of the mixture is given by 

l/g = P /g* A + (l~ p)/g* B + vp(l- P ), (4) 

where p is the number fraction of monomers A. Since the number of con- 
tacts per polymer is proportional to the density, the composition dependence 
of g translates directly into a composition dependence of the x parameter 
in blends of polymers with dissimilar monomer structure, g* A 7^ g* B . Thus 
X should depend linearly on the composition of the blend. The linear con- 
tribution vanishes in blends of monomers with very similar monomer struc- 
ture, e.g., isotopic blends or saturated hydrocarbon mixtures, and the re- 
maining composition dependence is parabolic^]. This is indeed observed 
experimentally jh^, |dj. Equation-of-state effects are thus apparently respon- 
sible for the linear part of the composition dependence of x- However, the 
magnitude of experimentally observed parabolic contributions cannot be ex- 
plained by Eqn.(|) alone 

A second important factor is the local structure of the fluid, i.e., the form of 
the correlation functions. Let us first note that the demixing is mainly driven 
by the miermolecular contacts between monomers of different chains. The 
miramolecular contacts contribute to the conformational free energy of single 
chains, which does not change very much upon mixing, since conformational 
changes are generally not very high]44|]. In the simplest approximation, the 
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pair distribution of monomers i and j from different chains is assumed to have 
the form 

Q { $ Mter {r) = QiQ j 9% ter {r) with gff "(r) = g inter (r), (5) 

where the pair correlation function g inter (r) depends neither on the type i 
and j of the monomers nor on the composition. It is normalised such that 
g(r) — ► 1 for r — > oo. In real fluids, this approximation will fail in two re- 
spects: First, the local packing of chains depends on the chain species, either 
directly due to monomer structure differences (monomer size etc.), or as a 
more subtle result of chain structure differences (chain architecture, chain 
stiffness). For example, the position of the peak in the correlation func- 
tion which corresponds to the first coordination shell depends on the size of 
the central monomer |^5j. The reasons for non-random packing are generally 
entropic: Excluded volume effects (related to the effective monomer sizes), 
effects of chain conformational entropyjl6| etc. Energetic interactions usu- 
ally do not affect the packing very much. A particularly strong effect on 
the demixing behaviour can be expected if like monomers pack closer than 
unlike monomers. Monte Carlo studies have shown that such "non-additive 
packing" alone is sufficient to bring about phase separation [47| . Second, 
even if the local fluid structure is preserved, i.e., the sum Y^j Qi 3 l '' mter (r) / gt is 
independent of i, one still expects local composition fluctuations. Such "non- 
random- mixing" also affects the demixing behaviour especially very close 
to the critical point Es]. According to the Ginzburg criterion [Ell, however, 
the random-mixing approximation becomes better upon increasing the chain 
lengthg9|. 

To summarize, demixing in polymer blends occurs for energetic and en- 
tropic reasons. The energetic factors include: Energetic incompatibility of 
monomers, and shifts of the ratio between inter- and intramolecular monomer 
contacts, caused by conformational changes. The entropic factors include: En- 
tropic incompatibility of monomers {e.g., due to non-additive packing), pack- 
ing inhomogeneities due to the different chain structure, and conformational 
changes of the chain. A huge amount of theoretical work has been devoted to 



elucidate the effect and the importance of the different contributions Eg, |31 
jspf . In most cases, the energetic or entropic incompatibility of monomers dom- 
inates the demixing behaviour. On the other hand, details of the monomer 
interactions are irrelevant on the scale of whole chains. It is thus reasonable 
to retain the spirit of the Flory theory and absorb the microscopic details 
into a few effective parameters, e.g., the x parameter and the compressibil- 
ity. These can be used as input parameters into theories of more complex 
systems such as polymer interfaces and surfaces. In addition, they allow to 
relate simulations of coarse-grained polymer models to experimental systems 
and idealised theories. 

The problem remains to calculate x f° r a given simulation model. Obvi- 
ously, Eqn.(H) cannot be used for any model different from the Flory lattice 
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model. However, the extension of the Flory theory to continuous space mod- 
els or more sophisticated lattice models is relatively straightforward. This 
shall be demonstrated in the following. We emphasise that we do not aim to 
present a complete theory of the \ parameter, nor to review the state of the 
art of generalised Flory theories. Rather, we wish to present a simple "recipe" 
for the calculation of \, one which takes into account the dominant contri- 
butions, and thus gives good results for most practical purposes. In fact, our 
type of approach has usually provided a good quantitative understanding of 
simulation data in the past jf6|, ^l], ^2|, |53| . 

We consider a mixture of ua, ns polymers of length Na, Nb- Polymers 
k are characterised by their center of mass position Rk and the relative coor- 
dinates of the monomers Uj : k = — R^. The general partition function of 
such a system in the given volume V can be written in the form 

(V/Vn) nA+nB 

Z= {V/V ° exp(-m/p,g)), (6) 

n A \n B \ 

where m = uaNa + ubNb is the total number of monomers in the system, 
Vq is an arbitrary reference volume, and f(p, g) is defined by 



nA+nB f, ft N, 



e -m f(p,g)) 



k=l Jn 7 = 1 



with (3 = 1/ksT and the total energy U. Note that the center of mass 
positions R^ have been rescaled such that the integration volume f£ does not 
depend on the volume V any more. The Helmholtz free energy thus reads 

0F = n A H^f) + n B In(^fr) + m f(p, g). (8) 
V/Vq v/Vo 

The first two terms describe the combinatorial entropy of mixing and the 
translational entropy of the center of mass of polymers. The last term sub- 
sumes the remaining contributions to the free energy, i.e., the internal energy 
and the conformational entropy. Since both are proportional to the number 
of monomers m, it is conveniently expressed in terms of a "monomer free en- 
ergy" f(p, q), which depends on the total density g and the number fraction 
of A monomers p. The pressure p at volume V is given by 

and the compressibility k by 

^-1 = ^ = 2^+^0, (10) 

where terms of order 1 /N (ideal gas contributions) have been neglected. Equa- 
tion (0) can be used as an implicit expression for the density g(p,p), as a 



8 



function of composition p and pressure p. At constant pressure, the Gibbs 
free energy per monomer p, then reads 

= tt m hr7-) + ln ( r T ) + AP> e) + — ■ n ) 

iVA iVA iVs iVs g 

One obtains the excess free energy of mixing 

/? M oxc = pn(j>,p)-ppn(l,p)-(l-p)l3ft(0,p) (12) 

= ^ln( /5 4) + i-^ln((l-p)4)+/3^ oxc 
N A q a n b Q b 

+[/(p, Q) - P /(I, A) -0--P) /(O, 0s)] (13) 

with the densities of the pure component melt f>l and g^, and the excess 
volume on mixing per monomer u cxc . M° xc can be used to define an effective 
Flory Huggins parameter x e ff 

/3M CXC = In ^ + ^ ln(l - 0) + Xeff P(l " P). ( 14 ) 
where <p is the volume fraction of monomers A. 

^4/(4 + ^)- d5) 

Qa 1 Qa Qb 

The temperature dependence and the phase behaviour results from the inter- 
play between the enthalpic and entropic contributions. Since the entropy of 
mixing is proportional to the inverse chain length, entropic contributions of 
different origin to the effective Flory-Huggins parameter \eS are much more 
important in macromolecular blends than in mixtures of small molecules. 

In order to proceed, it will be useful to expand f(l 7 g A ) and /(0, Q B ) 
around g. Using Eqns.(||) and (|To|), and dropping terms of order 1/N, this 
yields 

f(p, £?*) « f(p, Q) + ^(l - 4) + A(2/?P - /8«r X )(l - 4) 2 - (16) 

with the compressibilities of the pure component systems «j. Furthermore, 
we define the difference between volumes per monomer in the pure A and B 
melt, 8 = l/g*A ~ 1/ 1 Q*bi and assume that both 8 and the volume change upon 
mixing v cxc = 1/g— pg* A — (1 — p)/ 'g^ are small. To second order in gS and 
£w oxc , one obtains 

/3M exc = ^ In + ^ ln(l - 0) + A/(p, g) + /?/x^ p (17) 

with Af( P , g) = [/( P , e ) - p /(i, e) - (l - p) f(0, g)} 
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and the cquation-of-state contribution 



/z exc 



comp 



^gS 2 p(l-p) [(1 - P )k/ + P k b x ] 

+ \q v Ikc [PI^A 1 + (1 - 

+ £?W P (l - p) [Kg 1 ~K-/]. (18) 

So far, these are general thermodynamic considerations. The next step is 
to evaluate Af(p 7 g) and v exc . It is at this point that we introduce the gen- 
eralised Flory Huggins approximation. First, we adopt the Flory assumption 
(i) and neglect conformational changes upon mixing; hence intrachain inter- 
actions do not contribute to the free energy. Second, we assume that species 
dependent interactions between monomers can be treated as perturbations of 
a reference system "0" of polymers made of identical monomers [j54| . 

We suppose that we can choose a reference system with a short range 
monomer potential uo(r), e.g., a hard core potential, which does not dis- 
tinguish between monomers of different type. We further suppose that the 
correlation functions gl™ ter, °(r) for interchain contacts and gl" tra, ° (r) for in- 
trachain contacts, and the "bare" incompatibility Af°(p, g) in the reference 
system, are known. The correlation functions are related to the monomer pair 
distribution g^ 2 \r) via 

(2)/ \ 2 inter,0 / \ , c intra,0 1 \ 

Qij (r) = g PiPj g l3 (r) +gp i S ij g u (r), (19) 

where Sij is the Kronecker symbol. They arc independent of the composition 
p in symmetrical systems, but they may have a (weak) composition depen- 
dence if the polymers A and B have different architecture or stiffness. In 
practice, A/°(p, g) is very small and can usually be neglected, and gfj(r) can 
be approximated by the correlation functions gij (r) measured in simulations 



of the "full" system |§ [3§. 

Next, monomer specific "nonbonded interactions" are turned on. One has 
to distinguish between two possible factors: The monomers may have different 
size, i.e., the excluded volume interactions may depend on the type of the 
monomers, and additional energetic interactions may be present which also 
distinguish between monomers. The excluded volume part of the interaction 
will be strong, but short ranged, whereas the energetic part may be extended, 
but weak. It is thus advantageous to separate the total monomer potential 
Uij{r) into three parts 

Uij(r) = u (r) + Vij{r) + w i} {r), (20) 

where uq(t) is the potential of the reference system between identical monomers 
of radius d, Vij (r) is short ranged and nonzero only in a small region £d around 
r = d, and Wij (r) is weak. We will borrow concepts from perturbation ex- 
pansions in the theory of simple liquids ]5q| and proceed in two steps: First, 
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we include the potential Vij (r) and expand to lowest order in £. Then, we use 
this system as the new reference system and expand to lowest order in Wij . 

In order to do so, we introduce the indirect correlation function for inter- 
chain correlations [p5| 

y W(r)=^°(r)exp03uo(r)) ) (21) 

and, analogously, the indirect correlation function for intrachain correlations 
2/|"* ra ' (r). In simple liquids, Uijir) has the advantage of being continuous 
even in hard core systems. We shall assume that this holds also for polymer 
melts. In systems of polymers with equal stiffness, in particular, the interchain 
correlation functions in the polymer melt can be related to the corresponding 



monomer correlation function j/ m0 no('") via fl|, 31 > 4Sfl 

y% ter '°(r) = feono(r) (l - ifo (-^)) , (22) 



where gij are smooth functions (see also Sec. 4). If the monomers have hard 
core interactions, yij (r) outside the core is simply the pair correlation function, 
and good estimates for the value of y mon o inside can be looked up in the 
literature |55|. We now consider a "primed" intermediate system, in which 
Vij(r) is turned on, and uiij(r) is still set to zero. The perturbation expansion 
to lowest order in £ yields the monomer free energy 

f'(p, q) = /V Q) + \jdr {e-^W - e -^M+^«)} 

r \ ^ inter. / \ , \ A intra, 0/ \ /no\ 

[QZ^PiPjVij K r )+Z^PiVu M (23) 



\^ y _ nil' r.O, j \ r /iilrn. N, 

i,3 



and the indirect correlation function y[j{v) ~ Vijix) + C(£ 2 ) ( c f ©)■ Here 
the sum runs over i,j = A, i? and pi is the number fraction of component i, 
i.e., pa — p and pb = (1 — p). Corrections due to the chain connectivity J56| 
are of order 1/a/ZV and will be ignored here. The pair correlation function in 
the primed system thus takes the form g'^r) = yfj{r) exp[— f3(uo(r) + (r))]. 

We use the primed system as the new reference system and expand /(/?, g) 
in lowest order of Wij (r) 

f(p, Q) = f(p, 0)+^ /^«(r){e2piP i fl{/ nl ^(r)+2ft^ i totra (r)}.(24) 

i,j i 

After putting everything together, we get to lowest order 

f(p, Q) = f°( P , Q) + | £ ^ t£r (P) + \ E fX -Vu tra (P) (25) 

with effective interaction parameters jij (p) 

Hj(p)= / dry^(r){exp(-^uo(r) - exp(-/3u ii (r))}. (26) 
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The latter may depend weakly on the composition p via the pair correlation 
function (r) . It will prove useful to introduce the quantity 

X = 2 1 ^ er (p)~lAA r (p)~lBB r (p) (27) 

+ [{1AA\P) - 1AA 1 (1))/(1 -P) + {1BB 1 {P) ~ lBB l m/p\, 

where l\° tal {p) is defined as the sum i\° tal {p) = llf er (p) + l\f ra {p) I Q- Note 
that the bracketed term vanishes in symmetrical mixtures. In general, it is 
small and approaches a constant as p — > or p — > 1. 

We are now in the position to calculate the volume change upon mixing 
u xc and the volume difference S between the two pure component systems. 
Eqns.(p5|) and ( fL0| ) inserted into Eqn.(|j]) yield the composition dependence of 
the density 

q = eo (i - yb [*o E ftft 7™ ter (p) + eo I> t« H ) ' (28) 

where go is the density in the reference system. Using this expression, we 
obtain 

S=\-^ « ^(7^(0)-7il a '(l)) (29) 

v cxc = ; — ~ —Xp(l-p). (30) 

Q Qa Pb 2 P 

In particular, we recover the roughly parabolic composition dependence of v exc 
announced above. Note that these equations relate equation-of-state effects 
at constant pressure to quantities which are accessible in constant-volume 



simulations. Dudowicz and Freed 28 have studied the volume change of 



mixing in detail for a lattice model within lattice cluster theory. They find 
that fexc is slightly asymmetric, but the approximation by a parabolic law 
is still reasonable. Unfortunately, they do not calculate the compressibility 
and the correlation functions, thus we cannot relate Eqn.(pO|) quantitatively 
to their work. 

These results can be inserted into Eqn.(p^). The effective Flory Huggins 
parameter then takes the form 

Xeff = I * + + w [(7ir (0) " 7ir (1))2 + x2 p[l " p)m) 

The first term comprises the incompatibility of monomers and usually dom- 
inates by far. It contains equation-of-state effects in part via the prefactor 
g. The second term results from packing inhomogeneities due to the different 
chain structure, and the last term subsumes additional compressibility effects. 
Effects of conformational changes on mixing are disregarded, and so are non- 
random-mixing effects, which would be of higher order in the perturbation 
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expansion. Since the monomer contribution is generally much stronger than 
the other, the simplified version of Eqn. ( |3l"| ) 

Xeff * f X (32) 

yields a good estimate of the Flory Huggins parameter in most cases. 

Most simulations as well as analytical self-consistent field |2l], |2^, |2J| and 
P-RISM|3l| calculations are performed at constant volume. Here, the situa- 
tion is slightly different: The free energy of mixing is the difference between 
the Helmholtz free energy of a mixed system at density p, and between the 
sum of free energies of the corresponding pure component systems at densities 
g* A and g* B , where the average density is g, but the two pure systems are at 
pressure equilibrium with each other. In other words, the pressure of the pure 
systems p* is different from that of the mixed system, p, and chosen such that 
there is no volume change on mixing. Thus Eqn.@ has to be replaced by 

*~ = ^ln0 + I-^ln(W) (33) 
m Na Nb 

+ [f(p,g)- P f(l,Q A )-(l-p)f(Q,Q* B )}, 

where g* A = g(p = l,p*) and g* B = p(p = 0,p*) are evaluated from Eqn. ([)]), 
and the pressure p* is determined implicitly through the constraint 

4- + ^ = -- (34) 
Qa 8b Q 

The effective \ parameter can then be calculated in the same way as before (in 
Sec. 2.1). One obtains an expression similar to (|3l|), with different compress- 
ibility contributions. The situation is particularly simple for mixtures with 
8*A = 6*b- Eci n (HH) then implies g = g* A = g* B , and the effective x parameter 
is given by 

g Af°(p,g) 

^ = 2^ + 7M' (35) 

We will present some applications of these considerations to simulation models 
in section 4.1. 



3 Models and technical aspects 

Polymeric materials are characterised by widely spread length and time scales ]5j| 
[58| , |5{| . They range from the size of the chemical subunits (about A) to the 
coil size of the order 100A. In mixtures, the correlation length of composition 
fluctuations becomes even larger very close to the critical demixing point. Fi- 
nally, the length scale which characterises the morphology of an immiscible 
blend is macroscopic, and crucially determines the properties of the material. 
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The corresponding time scales range from the vibrations of chemical bonds 
(10~ 14 s) to the time needed for macroscopic phase separation, i.e., hours, 
days or even longer. 

A unified theoretical treatment spanning the whole range of length scales is 
clearly not feasible. Fortunately, polymeric materials share a variety of differ- 
ent dynamic and static properties, which can be traced back to a few relevant 
polymeric properties. This justifies the description of polymer mixtures on 
a coarse-grained length scale, using appropriate coarse-grained parameters. 
One of them is the x parameter, which has been discussed in the previous 
section, others will be introduced shortly. On that level, the study of coarse- 
grained polymer models can contribute to a better understanding of the blend 
properties in two ways: 

First, one needs to establish the relation between the chemical structure 
of polymers and the coarse-grained parameters. Unfortunately, a quantita- 
tive prediction of, e.g., the Flory-Huggins parameter x> from the atomistic 
structure is extremely difficult. This is because the typical energy scale for 
intermolecular interactions on the monomer scale is of the order /c^T, whereas 
the Flory-Huggins parameter is typically 3-5 orders of magnitude smaller in 
dense polymer blends. Here, coarse-grained models which incorporate some 
degree of fluid- like structure or conformational asymmetry provide useful qual- 
itative information about the relation between the microscopic structure and 
the coarse-grained parameters. 

Second, the models can be used to study the interplay between the coarse- 
grained parameters and the bulk thermodynamics and interfacial structure. 
Such investigations yield valuable insight into the universal material properties 
in polymeric systems. They can be compared to predictions of mean-field 
theories or experimental results. This aspect is the focus of most simulations 
presented here. 

A minimal set of relevant polymeric properties comprises: The connec- 
tivity of the macromolecules along the backbone, the excluded volume of 
the segments, and short range thermal interactions. Coarse-grained models, 
which retain these macromolecular properties, have proven extremely effi- 
cient in investigating the universal thermodynamic properties of polymeric 
multicomponent systems. Only a small number of parameters is required to 
compare experiments, simulations and theory quantitatively on this coarse- 
grained level. 

These are, for example, the \ parameter, the compressibility, and the chain 
length. In inhomogeneous polymer melts an addition parameter is required 
which sets the length scale in the system. A natural choice is the extension of 
single chains, i.e., the end to end distance R e of a bulk chain, or equivalently, 
the gyration radius R g . In a dense melt, sufficiently long and flexible chains 
are well described by random walks Q. Thus R e and R g depend on the degree 
of polymerisation N like R e = bVN and R g = b-i/iV/6, where the proportion- 
ality constant b is called the statistical segment length. The two parameters 
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xN and R e are generally sufficient to describe inhomogeneous systems of long, 
flexible polymers with weak monomer interactions (Gaussian chain model). 
Here, "weak" means that the length scale for composition changes is sev- 
eral times the statistical segment length, i.e., x "C 1. If the chains are very 
stiff [Q, or the incompatibility is quite large ]55|], a second length scale may 
come into play, the persistence length £, which characterises the length of the 



chain over which the monomers are still strongly correlated[100 



A minimal polymer model, which retains the above properties, is the lat- 
tice model of Flory and Huggins |26| . A small number of chemical subunits 
is represented by a single lattice site on a simple cubic lattice. The excluded 
volume property is mimicked by the constraint of single site occupancy, and 
bonded segments occupy nearest neighbour sites on the lattice. Phase sepa- 
ration is caused by unfavourable interactions between neighbour segments of 
different type. In spite of its extreme simplicity, this model and its mean-field 
analysis by Flory and Huggins is usually the basis for discussions of polymer 
miscibility in melts, and it already yields many important informations on the 
bulk thermodynamics (see section 2.). In particular, the Flory Huggins the- 
ory predicts a linear increase of the critical temperature of a binary polymer 
blend with the degree of polymerisation, in agreement with experiments J60| 
and simulations p8[|. 

Lattice models which resemble the Flory-Huggins-model have been used 
in simulations of polymer melts and solutions. One prominent example is the 
Larson model, which has found wide application for the study of amphiphilic 
systems (for a recent review seepi]). 

However, many important properties of polymer melts cannot be repre- 
sented in such a simple model ^nj: The small number of bond angles does 
not allow for a realistic modelling of the bending rigidity or other structural 
asymmetries. Monomers and vacancies have the same size, hence the un- 
derlying monomer fluid displays no packing effects [p2[ . These are however 
crucial ingredients to the Flory-Huggins parameter, as discussed above. With 
the constraint that the bond length is always one lattice unit, the implemen- 
tation of a local diffusive monomer dynamics is not possible. The latter is 
needed when studying the collective dynamics of phase separation (i.e., spin- 
odal decomposition), or the relation between the thermodynamic state and 
the single chain dynamics. 

The bond fluctuation model (BFM) of Carmesin and Kremer|33| over- 
comes the above mentioned limitations while retaining the computational ad- 
vantages of a lattice model: Each monomer blocks eight corners of a whole 
unit cell in a simple cubic lattice from further occupancy. Monomers along a 
polymer are connected by one of 108 bond vectors of length 2, \/5, y/6, 3, or 10. 
Due to the size difference between vacancies and monomers, one gets fluid-like 
packing effects. The extended monomer size and the variety of possible bonds 
helps washing out the underlying lattice structure. A local monomer hopping 
dynamics can be implemented, which leads to a Rouse-like dynamics in dilute 
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solutions. 

At the volume fraction 0.5 or monomer number density g — 0.5/8 (the 
factor 8 reflects the monomer volume), the model mimics many properties of 
a concentrated polymer solution or melt. The excluded volume interaction 
is screened down to about 6 lattice spacings, i.e., the chains have Gaussian 
statistics on larger length scales. The diffusion of very long chains exhibits 
the signature of reptation-like motion: the entanglement length is roughly 32 
monomers |j3. Due to the relatively high number of remaining vacancies, the 
polymer conformations still equilibrate reasonably fast. Real polymer systems 
can be mapped onto the BFMj65|, One monomer then corresponds to to 3- 
5 chemical repeat units, and one lattice units to roughly 2 A. In the following, 
all lengths shall be given in units of the lattice constant. 

The BFM has been widely used to study the statics and dynamics of 
polymer solutions f67|, |64| , melts f64|, |68| , and glasses f^. Hence, a variety of 
properties of the model are known, and this established model is well suited 
to test theories of polymer miscibility and polymer interfaces. 

Due to the high number of possible bond vectors, one has 87 different bond 
angles. This allows for a rather realistic modelling of a local stiffness along the 
chain. Unless noted otherwise, all of the polymers are complete flexible. In 
some of our studies [[l(| |53| , we impose on the B component of the blend an in- 
tramolecular potential, which favours straight bond angles: E(<3) — /cos(0), 
where O denotes the complementary angle between successive bonds. On in- 
creasing the strength of this potential from / = (flexible) to / = 2, the 
chain extension increases by a factor 1.5 (for chain length N = 32). 

We shall consider binary polymer mixtures containing Ha A-polymers of 
length Na and tib -B-polymers of length Nb, and ternary systems with ad- 
ditionally nc symmetric AB diblock copolymers of length Nc- In some ap- 
plications, the stiffness of the B component is tuned through the bond angle 
potential described above. The two kinds of monomer species, A and B, in- 
teract via sl square well potential. The range of the potential corresponds 
roughly to the first coordination shell of the monomer fluid, and includes 
the 54 nearest lattice sites. Unless noted otherwise, monomers of the same 
type attract each other, whereas contacts between unlike species increase the 
energy 

- e-AA = -£bb = cab = e. (36) 

In the following, we measure all energies in units of ksT. Hence, e is propor- 
tional to the inverse temperature. 

One might argue that purely attractive monomer interactions of different 
strength would mimic the experimental situation in a more realistic way. The 
choice ( |36|) is motivated by the following consideration: Even the pure one 
component system (p = or p = 1) is in fact a mixture of homopolymers and 
vacancies. For strong attractive interactions between monomers, one hence 
encounters a liquid-vapour phase separation into a homopolymer-rich and 
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vacancy-rich phase. The temperature of this "0 transition" is independent 
of the chain length (ee ~ 0.495(5) (67|), whereas the temperature of the AB 
unmixing transition increases like N. For long chains, these two tempera- 
ture scales are thus well separated. For short chains, however, the choice of 
purely attractive monomer interactions would lead to a rather complicated 
interplay between both transitions. It is thus expedient to choose the inter- 
actions such that the demixing occurs already at rather weak absolute values 
of eij. Simulations of the Bond fluctuation model have also been performed 
using a slightly smaller interaction range. The results for strictly symmetric 
interactions are qualitatively similar (||, |49| . 

Even though complex lattice models capture many of the relevant proper- 
ties of polymer mixtures, there are limitations: Constant pressure simulations 
of lattice models [fro] are difficult, whereas their implementation is relatively 
straightforward for off-lattice models. Keeping the pressure constant is de- 
sirable, for example, when the blend exhibits a pronounced volume change 
upon mixing (cf. section 2.). Even at constant volume, off-lattice models 
have the advantage that the pressure and the interfacial tension can be mea- 
sured via the virial equation fl47[ [7l| . Furthermore, they lend themselves to 
molecular dynamics simulations, which are able to capture hydrodynamic 
flow. The latter is crucial, for example, to describe the final stage of spin- 
odal decomposition]^]. On the other hand, off-lattice simulations are more 
demanding with respect to computational resources. In many respects, the 
BFM simulations are remarkably consistent with off-lattice predictions, e.g., 
with respect to the packing contributions to the \ parameter (see Sec.4.1). 

In studies of the thermodynamics and structure of polymer blends, two 
statistical ensembles are mainly used: In the canonical ensemble, the temper- 
ature, the volume and the global composition of the mixture are fixed, and 
the Monte Carlo moves update the polymer conformations. In simulations 
of the BFM, local monomer hops(6^] and slithering snake-like movements [fz3| 
have been applied. The former allow an interpretation of the Monte Carlo 
simulation in terms of an overdamped, purely diffusive dynamics, the latter 
relax the polymer conformations faster by a factor N. In off-lattice models, 
MD integration^] is often employed to propagate the system. The canon- 
ical ensemble is useful for the study of dynamic properties, the local inter- 
facial structure and the phase behaviour of systems with strong structural 
asymmetries [^3) . 

In systems with weak or moderate structural asymmetry, however, the 
phase behaviour is investigated more efficiently in the semi-grandcanonical 
ensemble [fZ4|| : The temperature, the volume, and the total number of monomers 
are held constant, while the composition of the mixture is allowed to fluctu- 
ate. In binary mixtures, it is controlled by just one parameter, the exchange 
monomer potential A/i (see Eqn. (fl9|). In ternary mixtures, a second variable 
S/i has to be introduced, which couples to the copolymer content in the mix- 
ture. The Monte Carlo scheme includes moves which turn polymers of type 
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A into polymers of type B and vice versa, and additional moves in ternary 
systems which switch between homopolymers and copolymers |7j| . This al- 
gorithm has first been introduced by Sariban and Binder j74j, and applied to 
the simple Flory-Huggins lattice model of strictly symmetric binary blends. 
In that case, the semi-grandcanonical identity switches are merely a change 
of the polymer label. It literally carries over to off-lattice models of sym- 
metric mixtures (47]]. The scheme can also be applied to asymmetric mixtures 
with different dispersion forces [{76) , with chain length asymmetry^, f73f , 
with stiffness disparity and non-additive monomer shapes f4(| . For even 
stronger asymmetries, configurational bias |67|| or gradual inserting schemes 
can be envisaged. For very strong conformational asymmetries, however, the 
semi-grandcanonical moves become increasingly inefficient [ |53| . 

The advantages of the semi-grandcanonical ensemble are threefold f75f : 

First, since the global composition is not conserved, large length scale com- 
position fluctuations relax much faster than in the canonical ensemble. The 
attempt to observe an unmixing transition driven just by the diffusion of poly- 
mers poses a substantial challenge in terms of computational resources ]6l[ . 
Furthermore, a thorough analysis of finite-size effects in the canonical en- 
semble {e.g., via a subbox analysis ]77[) requires formidably large system 
sizes. The semi-grandcanonical ensemble concentrates the computational ef- 
fort on the composition fluctuations, the slowest mode in the blend. Sophisti- 
cated finite-size scaling methods have been applied to study the crossover 
from mean-field to Ising critical behaviour in strictly symmetric polymer 
mixtures [Q, and mixed-field finite-size techniques have been used to accu- 
rately determine the critical point in asymmetric mixtures p7|, [78], f79j. These 
techniques can be combined with a multi-histogram extrapolation of the joint 
energy-composition distribution |6l], |8TJ and hence allow to determine the 
phase behaviour accurately from simulations of systems of rather modest size 
(i.e., comprising only a few hundred polymers). 

Second, the semi-grandcanonical equation of state, which relates the com- 
position p to the chemical potential difference Ap, can be measured in semi- 
grandcanonical simulations. This gives information about the thermodynam- 
ics of mixing, and provides the most direct connection to mean-field theories. 
In particular, it can be used to relate simulations to the considerations of Sec. 
2 (see also Sec. 4.1). 

Third, the semi-grandcanonical ensemble can be combined with multi- 
canonical reweighting schemes [jsT|| . Rather than generating configurations 
following the Boltzmann distribution, the simulation then samples the phase 
space according to the probability exp(— (E — 7iAAp)/kBT)/W(p), where E 
denotes the total energy of the binary system. By choosing the preweight- 
ing factors W proportional to the equilibrium probability distribution in the 
semi-grandcanonical ensemble P{p), one achieves uniform sampling of all com- 
positions p. An excellent estimate for the equilibrium probability distribution 
can be obtained from histogram extrapolations of previous runs. Below the 
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critical temperature, this Monte Carlo scheme|5l], 75 encourages the system 
to explore configurations in which the two phases coexist in the simulation 
box. 

We shall illustrate the last point with the simulation data for a symmetric 
ternary homopolymer/copolymer blend. The probability distribution of the 
composition as a function of the chemical potential S/i of the copolymers is 
presented in Fig{j] on a logarithmic scale. We have chosen an elongated simu- 
lation cell|7f| of xy-cross section L 2 and extension in ^-direction 2L = 96. The 
distributions exhibit two pronounced peaks corresponding to the composition 
of the two coexisting phases. In the two-domain state between these coex- 
istence compositions, the probability is strongly reduced. For large enough 
system sizes, the probability distribution there is dominated by configurations 
containing two interfaces parallel to the xy plane, as sketched schematically 
in the figure (periodic boundary conditions are applied). The two interfaces 
of area L 2 separate the two bulk coexisting phases. The probability of these 
configurations with respect to the bulk is reduced by the Boltzmann weight 
exp(— 2L 2 a /ksT) of the interfacial free energy cost|82 , where a is the inter- 



facial tension between the coexisting phases. Hence, the interfacial tension 
can be calculated from 

a 1 Inj^j) I--*") 



k B T 2L 2 V Pi 1 / 2 ) 

The plateau in the probability distribution around p — 1/2 shows, that the 
two interfaces can vary their mutual distance, thereby changing the compo- 
sition, without free energy costs. This important consistency check indicates 
that the the interfaces can be described as flat and non-interacting. The free 
energy of the interface is lowered upon increasing the copolymer chemical 
potential. Within the same simulation run, we can determine how the num- 
ber of copolymers at fixed chemical potential Sfj, depends on the composition, 
(cf. Fig. | (b)). Not surprisingly, configurations with two interfaces contain 
more copolymers. From this we can extract the excess of copolymers at the 
interface. It grows upon increasing the chemical potential Sfi of the copoly- 
mers. Similarly, we can monitor the energy as a function of the composition, 
and extract the excess energy of the interface. In addition, a careful anal- 
ysis of the probability distribution yields information about the interaction 
potential (joining pressure) between two interfaces. Hence, bulk and excess 
interfacial properties are simultaneously accessible in the semi-grandcanonical 
ensemble J75[. 

The same combination of techniques has been applied to thin film, where 
both flat walls attract the same component S3]. In this case the coexistence 
chemical potential A/i coex is shifted from its bulk value A/i^ oea . = 0, but 
nevertheless reweighting techniques permit to locate the phase coexistence 
accurately even for low temperatures. 

Even with a careful choice of the statistical ensemble and the use of recently 
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Figure 1: (a): Probability distribution of the composition in a ternary 
homopolymer-copolymer mixture at e = 0.054 and system size 48 x 48 x 
96. Upon increasing the chemical potential Sfi of the copolymer, the 
"valley" in the probability distribution becomes shallower, indicating 
that the copolymers decrease the interfacial tension, (b): Average num- 
ber of copolymers nc{p) as a function of the composition for the same 
system as in (a). 
From Ref. [fr5[ . 
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developed simulation and analysis techniques, the accurate measurement of 
structural and thermodynamical properties in the computer simulation poses 
huge computational demands. One the one hand, very large system sizes 
are required to investigate, e.g., the effect of capillary- wave broadening. In 



our simulations pi p6|, 117] we vary the lateral extension of the simulation 



box from 64 to 512 lattice units. The largest simulation cell contains more 
than a million monomers. Only massively parallel computers like the CRAY 
T3D/T3E make the simulations of such large systems feasible. Using a two- 
dimensional geometric parallelisation scheme, we have efficiently implemented 
the local monomer movements. Even for small simulation cells, the program 
scales practically linear with the number of processors up to 256 Q. On the 
other hand, many systems are characterised by protracted long time scales: 
Increasing the chain length N increases dramatically the relaxation times of 
the chain conformations. Deutsch et al. have employed chains of up to 512 
segments to investigate the crossover between mean-field and Ising critical 
behaviour. Their investigation required several thousand hours of CRAY Y- 
MP timeQ]. The equilibration of long wave length interfacial fluctuations is 
also difficult because of its slow dynamics. Another example is the relaxation 
time in the framework of the multicanonical reweighting scheme |8l[. In the 
most ideal case, the system performs a random walk on all compositions. 
Therefore the correlation time scales with the system size L like L 6 . The total 
amount of CPU time investigated into the results presented in the following 
exceeds the equivalent of 10 5 hours on a single CRAY T3D processor. 

In blends with strong structural asymmetries (e.g., large stiffness disparity, 
polymer solvent systems), the semi-grandcanonical identity switches between 
different polymer types are rather inefficient fl6|. [53[ . Hence, methods for mea- 
suring the interfacial tension in the canonical ensemble have been explored. In 
off-lattice models, the interfacial tension can be measured via the anisotropy 
of the pressure tensor. This method has been successfully applied to de- 
termine the free energy costs of a hard walljnj in a concentrated polymer 
solution and the interfacial tension |47j| in a binary polymer blend. However, 
the generalisation to lattice models is difficult ]7Q] . 

An alternative way of measuring the interfacial tension is the analysis of 
the capillary fluctuation spectrum [[53[ fn\ , [ssf : In general, interfaces are not 
flat, but exhibit long-wavelength capillary-wave fluctuations. A typical snap- 
shot of the local interfacial position in a binary polymer blend is presented 
in Fig.||. The local interfacial position has been laterally averaged over a 
size comparable to the polymer's radius of gyration, and large length scale 
fluctuations are clearly visible|^6|. Let us consider a flat interface, neglect 
bubbles and overhangs, and describe the local deviations of the interfacial 
position from its lateral average by the single valued function u(x, y). Due to 
the local deviations, the interfacial area is increased compared to the corre- 
sponding planar interface. The effect on the thermodynamics of the interface 
is commonly described by the capillary-wave Hamiltonian for long-wavelength 
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Figure 2: Typical snapshot of the local interfacial position u(x, y) for a 
system of lateral extension L — 64, e = 0.03 and N = 32. The coarse- 
graining size B = 8 corresponds roughly to the chain's radius of gyration. 
Long wavelength fluctuations are clearly visible. 
From Ref. p7|. 
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interfacial fluctuations |87|, Q 

H = J dxdy ||(Vu) 2 }, (38) 

where a denotes the effective interfacial tension. It can be diagonalised very 
naturally with a simple Fourier transform 

L/2-1 

U ( X )~Y + 51 { a (*) cos (<7z20 + Kit) sin (ll x )} > ( 39 ) 

1=0 

with qi = 2-kI/L. The equipartition theorem shows that the Fourier compo- 
nents have a Gaussian distribution with Gaussian width 

(40) 
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Figure 3: Probability distribution of Fourier components a(q) and 6(g) 
of the local interfacial position for the four smallest wave- vectors q in a 
system with system of lateral size L = 64, bending energies / = 0, 2, and 
chain length N — 32. For the data collapse the 1/q 2 dependence of the 
variance has been exploited. The solid line shows the expected Gaussian 
distribution. 
From Ref. Q. 

The probability distribution for the Fourier components of the local interfa- 
cial position are presented in Fig.0 for a symmetric blend and a blend with 
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strong stiffness disparity [£3|. In both cases, the long-wavelength interfacial 
fluctuations are well described by the quadratic interfacial Hamiltonian. The 
interfacial tension can be determined from the inverse width of the distribu- 
tion. For the symmetric blend, this value can be compared with the indepen- 
dent measurement in the semi-grandcanonical ensemble. The agreement is 
very good. The interfacial fluctuation spectrum thus provides an efficient al- 
ternative way of measuring the interfacial tension in structurally asymmetric 
systems. For example, the method has been successfully applied to measure 



the interfacial tension of amphiphilic bilayers|85 . 

We note that the relative increase of the interfacial area is proportional to 
h\L/L 2 . Hence interfacial excess quantities [f75| (e.g., the interfacial tension, 
the excess energy of the interface, the enrichment of copolymers or vacancies) 
are rather insensitive to capillary fluctuations. Interfacial profiles are affected 
more dramatically Q. "Apparent" profiles obtained via lateral averaging de- 
pend sensitively on the extension of the system and are significantly broader 
than the "intrinsic" profiles. This capillary-wave broadening has to be ac- 
counted for when comparing simulation results analytical predictions, and is 
also relevant for the interpretation of experiments. The fluctuations of the lo- 
cal interfacial position (u 2 (x,y)) are obtained by integrating over the lateral 
Fourier components. The power law spectrum in Eqn.(f40|) leads to logarith- 
mic divergencies for q — > and q — > oo. We can remove these divergencies by 
introducing heuristic cut-offs g max and q m im and obtain 

<" 2 (^)> = fe?lnf^V (41) 



27TCT \ q, 

The lower cut-off q m - ln is simply determined by the system size, i.e., q min = 
The value of the upper cutoff g max , however, is less obvious. In addition to 
long-wavelength capillary fluctuations, one has composition fluctuations on 
short length scales at the interface (i.e., intrinsic fluctuations), which deter- 
mine the shape of the intrinsic interfacial profile. The upper cutoff g max is 
hence related to the lateral length scale which separates the intrinsic fluctu- 
ations from the fluctuations of the center of the intrinsic profile. For small 
molecules and far from the critical region, one usually estimates <7 max = 27r/a, 
where a denotes the molecular extension. In polymeric systems pft, however, 
three possible microscopic length scales have to be considered: the statistical 
segment length b, the width of the intrinsic interfacial profile Wq, which is 
controlled by the Flory-Huggins parameter, and the radius of gyration, which 
scales like >/~N. Semenov |p9[ proposed to use g max = 2/wq as the upper 
cutoff. In principle, Monte Carlo simulations could address this problem by 
studying the dependence of the broadening on the chain length N and the 
incompatibility x- 

A finite bending rigidity n of the interface also generates a smooth cutoff of 
the fluctuation spectrum on small length scales. This effect can be described 
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by the Helfrich Hamiltonianj88| 

H = Jdxdy {|(Vu) 2 + |(A M ) 2 }. 

Including the bending rigidity k removes the divergency for q 
analysis as above yields 
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Thus the crossover length \/<j/k acts as upper cutoff q mayi . This effect is 
important in interfaces with absorbed amphiphilic molecules |75[. At high 
concentrations, they oppose themselves to being squeezed together on one 
side of a bent interface, thus generating a positive contribution to the bending 
rigidity (cf. Sec. 4.3). 

The laterally averaged "apparent" profile can be approximated by the 
convolution of the intrinsic profile po{z) with the Gaussian distribution P{u) 
of the local interfacial position 



p{z) = / du po{z — u)P{u). 



(44) 



For the slope at the center of the profile, this convolution approximation yields 
in the case of weak broadening pel 
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The slope at the center [z — 0) of the profile can be used to define the inverse 
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where we have approximated the shape of the intrinsic profile by an error 
function in the last step in order to calculate the numerical prefactor. 

Since analytic mean-field theories generally yield intrinsic profiles, one has 
to take due account of the capillary-wave broadening when comparing them 
with simulation data. In order to reduce the broadening effect, we can define 
"reduced" profiles by laterally dividing the system into subsystems of size 
B x s|5^, |86). In each subsystem, we localise the interfacial position and 
average profiles with respect to the local midpoint of the profile. In principle, 
the lateral length scale B can be chosen as to match the width of the reduced 
profile with the theoretically predicted value. However, once B is chosen in 
this way, other profiles {e.g., segmental orientations, chain ends, etc.) can be 
compared to theoretical predictions without any adjustable parameter. 
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4 Simulation results 



Many of the simulation and analysis techniques presented in the previous 
sections can be applied to lattice models of polymers as well as to off-lattice 
models, and also to non-polymeric models. In this section, we will show how 
they have been applied specifically to simulations of the bond fluctuation 
model, and compare the results to experiments and to analytical predictions. 



4.1 Local packing and miscibility 

We will begin with a discussion of the bulk thermodynamics of polymer mix- 
tures. The knowledge of the phase behaviour is required for a discussion of 
interfacial properties and of considerable interest in itself p8|. 

In the past, special emphasis has been given to simulations of symmetric, 
binary polymer blends p8[|. Experimentally, nearly symmetric blends can be 
realised by mixing partially deuterated and protonated polymers |}(|, or as 
blends of statistical copolymers Simulations by Deutsch and Binder p8| 
have confirmed the linear scalingp^] of the critical temperature by investi- 
gating polymers with up to 512 monomers. Using a sophisticated crossover 
finite-size scaling technique, they accurately located the critical temperature 
and investigated the crossover between mean-field behaviour to 3D Ising crit- 
ical behaviour, which prevails in the vicinity of the critical point. Upon in- 
creasing the chain length, the Ising critical region shrinks, in accord with the 
Ginzburg criterion |l8| and with experimental observations |60|, |92| . 

The phase behaviour of the bulk is determined by the local structure of 
the polymeric fluid, as reflected by the fact that the Flory-Huggins parameter 



X can be related to the intermolecular pair correlation function) 49 (see Sec. 
2.). In the simplest case of an athermal melt (i.e., e = 0) the intermolecular 
packing is driven by two effects: On the one hand, the size disparity between 
vacancies and monomers gives rise to fluid-like packing effects which result 
in a highly structured correlation function at short distances. On the other 
hand, the extended structure of the macromolecules causes a reduction of 
contacts with other chains on intermediate length scales. In order to separate 
the monomeric, local packing effects from the universal behaviour of the poly- 
meric correlation hole, we consider the ratio between the intermolecular pair 
correlation function and its monomeric equivalent. As shown in Fig. (^), it 
is found to be largely independent from the packing of the monomers. Since 
the correlation hole comprises N — 1 monomers and the only length scale of 
a flexible Gaussian polymer is its radius of gyration R g ~ y/N, the reduced 
correlation function exhibits the following scaling behaviour 

rr = 1 — y \ i — 

9N=x(r) r VviV, 

Fig.|] demonstrates for various chain lengths that such a scaling works 
very well for flexible chains. We note that the situation is somewhat more 



9 \-f= • (47) 
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complicated in systems of semifiexible chains, since the correlation hole is 
characterised by two length scales, the persistence length and the radius of 
gyration |||. 




Figure 4: Scaling of the reduced correlation function gjv( r ) /qn=i (r) with 
chain length. N = 20 (squares), N = 40 (diamonds), TV = 80 (triangles), 
and N = 160 (stars) for athermal (e = 0) melts. 
From Ref. §§. 

The deepening of the correlation hole with increasing chain length induces 
a chain length dependence of the effective number of intermolecular contacts 

~ C011St (48) 



N 

and hence presents a correction of relative order 1/yN to the scaling of the 
critical temperature for flexible chains (see below). 

In the following we will present examples for different contributions to the 
effective ^-parameter. In particular, we will consider the "semi-grandcanonical 
equation of state", which relates the composition of the blend p to the ex- 
change chemical potential A/i, 

/3A M = = — ln0 - — ln(l - <t>) - X e «W ~1) + C, (49) 

where C is a constant. Since this thermodynamic relation is directly accessible 
in the simulations of the semi-grandcanonical ensemble, we can extract an 



27 



effective x parameter. This value can be then compared to Eqn.(p5|), which 
relates the local fluid structure to the incompatibility x- 

Our first example is a blend of polymers of different chain length (Na 7^ 
Nb), which demix for energetic reasons, i.e., unlike monomers repell each 
other according to Eqn. (|36|) . Here the athermal (e = 0) system of non- 
interacting monomers suggests itself as the reference system. The bare contri- 
bution A/ to the x parameter is negligible, and the pair correlation functions 
in the reference system do not depend on the identity of the monomers. The 
combination of Eqns.(p6|), (f27|), and ( |35| ) then yields the effective Flory Hug- 
gins parameter 

Xeff = 2z c e, where z c = g dr g(r) (50) 

Jo 

defines an "effective coordination number" z c in the spirit of the original 
Flory-Huggins theory. We recall the g(r) is the intcrmolecular pair correla- 
tion function. FigJ^ demonstrates that the equation of state obtained with 
this value of Xeff ^ s m excellent agreement with the simulation data [fill , for 
temperatures well above the unmixing transition. 

The second example illustrates demixing driven by entropic repulsion be- 
tween monomers. One particularly simple realisation is non-additive packing. 
In real polymer mixtures, it might arise from disparities in the monomer shape 
or steric hindrances ]28|, . It is also used by Grest and coworkers to induce 
phase separation in an off-lattice model j|7]]. In our case, it is modelled as 
follows: Monomers of different type are subject to the additional constraint 
that they must not come closer than v5 lattice constants, whereas monomers 
of the same type can approach each other up to 2 lattice constants, as before. 
The natural reference system is again the athermal system of non-overlapping 
monomers. Eqns (|26|), ( p7| ) and (^) lead to the entropic x parameter 

Xeff = z & with Ze = Q drg(r). (51) 

Jo 

This non-additive packing thus gives rise to a rather large, positive contri- 
bution to the effective Flory-Huggins parameter, and the system demixes. 
In addition, energetic interactions can be introduced as before. The total 
effective x parameter is then given by 

Xeff = z 6 + 2ez c . (52) 

The case of attractive interaction between unlike monomers (i.e., e < 0) is 
particularly interesting. Here one encounters a Lower Critical Solution Point 
(LCSP) upon heating the system, i.e., reducing the strength of the thermal 
interaction e. 

Using the intermolecular pair correlation function of a reference system 
(|e| = 0.05, additive), we obtain z 6 = 0.238(2) and z c = 1.41(2) at the 
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Figure 5: Equation of state for an asymmetric polymer mixture Na = 10 
and Nb = 20 for various temperatures above the critical point, e = 0.01 
(circles), e = 0.02 (squares), e = 0.025 (diamonds). e c = 0.0320(1). The 
straight line corresponds to the prediction with z c — 2.735. 
From Ref. g| 
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monomer volume fraction 8g = 0.35, which corresponds to a concentrated 
solution. FigJ| compares the effective \ parameter extracted from the semi- 
grandcanonical equation of state in the one phase region {i.e., at low tem- 
peratures) with the above estimate. Eqn.(E)2|) describes the simulation data 
quite successfully ]46[ ]. We note that Eqn. ( j52p predicts that the critical tem- 
perature of the LCSP does not increase upon increasing the chain length, but 
converges to 1/T C = |e c | = z e /2z c . This is indeed found in the simulations, 
and in qualitative agreement with experiments of Schwahn and co-workers |94| 
on the phase behaviour of PVME/PS blends. 



0.2 



-0.4 



-0.6 - 



-0.8 







i 


i i 1 i 

0.2376-2.82|e| 
























O 






i 


o 

1,1,1 



0.00 0.05 0.10 0.15 0.20 0.25 
lei 



Figure 6: Effective Flory-Huggins parameter extracted from the semi- 
grandcanonical equation of state for a non-additive polymer mixture of 
chain length N = 20 at a total monomer density of 8g = 0.35. There is an 
attraction between unlike monomers (e < 0). Upon heating i.e., reducing 
the thermal interactions, the blend phase separates. The location of the 
Lower Critical Solution point is indicated as a vertical line. 
From Ref. @. 

The last example presents a situation where polymer incompatibility is 
caused on the level of whole chains as a result of different chain stiffness. 
Monomers are taken to be identical otherwise. Hence we are left with the 
"bare" contribution to the x parameter, \ — A/°(p, (?)/[p(l — which 
results from packing inhomogeneities due to the different chain structure. The 
explicit evaluation of this term requires more elaborate techniques |37], |3^] than 
introduced in Sec. 2. Fig.0 shows the semi-grandcanonical equation of state 
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for an athermal (i.e., e = 0) mixture of flexible (/ = 0) and semiflexible (/ = 
1) polymers. A composition dependence in the figure indicates a non-zero 
Flory-Huggins parameter. The simulation data clearly show, that the purely 
entropic packing differences give rise to a small, chain- length independent, 
positive Flory-Huggins parameter x- Though the effect is quite small, it 
might lead to phase separation for very long chain lengths N > 0(1000). 
Of course, one could suspect, that these subtle packing effects are strongly 
influence by the structure of the underlying lattice. However, field theoretical 
calculations by Liu and Fredrickson[[58) as well as recent P-RISM calculations 
by Singh and Schweizer[[37| predict effects of similar magnitude within off- 
lattice models. Off-lattice simulations by Weinhold et al. |S also confirm the 
order of magnitude of the packing induced incompatibility. 
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Figure 7: Semi-grandcanonical equation of state for a polymer mixture 
of flexible (/ = 0) and semi-flexible (/ = 1) polymers without thermal 
interactions (e = 0). The translational entropy is already accounted for 
and the composition dependence indicates the stiffness-induced contri- 
bution to the Flory-Huggins parameter. 
FromRef. M. 



Since purely athermal systems of flexible and stiff chains with chain lengths 
accessible to simulations usually don't demix, additional, repulsive monomer 
interactions have to be turned on in order to induce demixing. The "ather- 
mal" system of flexible and stiff chains then has to be used as the reference 
system |5i| (cf. Eqn.(p6|)). We note that the reference system is asymmet- 
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ric, and the effective coordination numbers depend on the identity of the 
monomers and the composition of the blend. This is because stiffer chains 
are more extended, and folding back is less probable. Hence the number of 
intermolecular contacts increases and, as a result, the interaction parameter 
7i, in ( p6| ) is larger for the stiffer species, even if the interaction potentials are 
chosen symmetric. The effect on the x parameter is comparable in magni- 
tude to the stiffness induced (athermal) contribution discussed above Q for 
temperatures of the order T c . For strongly incompatible blends, the enthalpic 
contribution to x dominates. 

Fig.|| shows an interface between two demixed phases in a blend of flexible 
A and stiff B chains, confined between hard walls. The two phases are at 
pressure equilibrium, and vacancies are enriched in the stiffer phase. A similar 
effect has been observed in an athermal two-dimensional systemj95|. It is due 
to the fact, that the stiffer chains pack less efficiently and hence their osmotic 
pressure is slightly higher. The order of magnitude of the effect is however 
only 1 %. The oscillations of the density directly at the walls are packing 
effects and reflect directly the monomer correlation functions. 
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Figure 8: Polymer blend of flexible (/ = 0) and semi-flexible (/ = 1.5) 
polymers (N — 16) between two hard and neutral walls well below the 
critical temperature T — 0.346T C - The normalised monomer density 
profile shows pronounced packing effects at the wall. The blow-up shows 
that vacancies are enriched in the stiffer phase (right hand side). 
FromRef. §§. 
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Simulation studies of polymer compatibility in continuous space and at 
constant pressure have been performed recently by Kumar |9(|. He studies 
miscible and immiscible blends of Lennard-Jones chains, and finds volume 
changes of mixing with the expected parabolic composition dependence. He 
also measures incremental chemical potentials of both components and re- 
lates them to an expression similar in spirit to Eqn.(32). Inserting correlation 
functions which include both interchain and intrachain contacts gives worse 
results than using correlation functions for interchain contacts only, in agree- 
ment with the considerations of Sec. 2. In particular, he studies the special 
case of a blend with interactions chosen such that X happens to cancel to 
zeroQ. As expected, the volume change of mixing practically disappears. 
Nevertheless, the blend demixes, if the chains are sufficiently long. The effect 
can be attributed in part to compressibility effects (i.e., the remaining third 
term in Eqn. (|3l|)), but probably also involves higher-order effects, such as 
non-random- mixing and chain-stretching effects. A systematic investigation 
of the scaling of the critical temperature with chain length would help to 
separate the different contributions, but has not yet been attempted. 

In the above examples we have related the structure of the polymeric 
fluid to its thermodynamical properties far from the unmixing transition, and 
obtained a reasonable quantitative agreement. Within mean-field theory, the 
unmixing temperature is given by Eqn.^|. 

However, even in the simplest possible systems - strictly symmetric, binary 
polymer blends - mean-field theory overestimates the critical temperature at 
short chain lengths [|49|. This discrepancy between the mean- field theory and 
the Monte Carlo results can be attributed to composition fluctuations. The 
fact that we observe 3D Ising critical behaviour at the unmixing transition 
already indicates that the deviations from the random-mixing approximation 
involved in the mean-field theory are quite strong. The effect is well known 
in mixtures of small molecules, i.e., the small chain length limit of polymer 
blends. P-RISM calculations |5fJ and calculations of fluctuation corrections 
to mean-field by Holyst and Vilgis(98| suggest that these non-random-mixing 
effects die out with l/v^/V- Thus, the ratio between the actual critical temper- 
ature and the mean- field value should converge to 1 with a 1/ VN correction. 
This is illustrated in Fig|| for strictly symmetric models with different inter- 
action range, and for blends with chain length asymmetry [fl9| . The relative 
deviation between the critical temperature and its mean-field estimate is the 
larger, the smaller the interaction range. When the mean-field theory is used 
to fit Monte Carlo data on the critical isotherm, non-random mixing effects 
thus lead to a strong parabolic composition dependence of the effective x 
parameter close to the critical point Es, Eg]. 
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Figure 9: Chain length dependence of the ratio between the critical 
temperature of a binary polymer blend and the mean-field estimate, 
which takes the local fluid structure into account. Symmetric mix- 
tures (Na = Nb) are presented by circles and triangles (for a re- 
duced range of the thermal interactions). Asymmetric mixtures with 
Nb — kNA are shown as squares (k = 2) and diamonds (k = 3). 

A/red = 4N A N B /(N A /2 + N'J 2 ) 2 . 

From Ref. 



34 



4.2 Homopolymer interfaces 

We shall now discuss some general concepts on the structure of interfaces 
between demixcd homopolymer phases. For the sake of simplicity, we will 
consider a symmetrical system of polymers A and B with equal chain length 
./V and statistical segment length b. Generalised expressions for asymmetric 
systems will be given later if necessary. Within the Gaussian chain model, 
the behaviour of the interface is basically driven by the relative polymer in- 
compatibility X-/V. One distinguishes the two limiting cases xN ~ * 00 (strong 
segregation limit) and \N — > (x-/V) c (weak segregation limit). The two quan- 
tities in which we will be interested are the interfacial width and the interfacial 
tension. 

The intrinsic interfacial width w results from an interplay between entropic 
factors and "energetic" factors related to the monomer incompatibility, i.e., 
the x parameter. The interfacial tension, a, is closely related. It turns out 
to be roughly proportional to m 2 /w, where m — p c £ ex — p c £ ex is the width 
of the miscibility gap. In the limit of very long chains, the conformations 
of chains are strongly perturbed in the vicinity of the interface. Polymers 
can win conformational entropy by allowing chain portions to loop into the 
unfavourable side. The entropy gain due to the formation of such loops has to 
be balanced with the energy loss due to the unfavourable contacts of the loops. 
Since the total chain length does not enter this argument, the interfacial width 
is independent of chain length, and so is the interfacial tension. The strong 
segregation theory (2/N <C x 1) yields plf 



u>ssL = fr/v 6 X, fassh = V x/6 q b. (53) 

In the limit of short chains or weak incompatibility, the conformations of 
chains are only weakly perturbed. Here, polymers as a whole may move over 
to the unfavourable side, thereby gaining translational entropy. The latter 
has again to be balanced with the energy loss, and the resulting interfacial 
width and interfacial tension depend on the distance of %iV from the critical 
point (xN) c . In the weak segregation limit, one thus gets 

w>wsl oc (1 - XclxY v -, cwsl oc (1 - Xc/xY, ( 54 ) 

with the critical exponents v = 1/2, /i = 3/2 and the critical point Xc = 2/iV 
in the mean-field regime. Experiments are usually well described by mean-field 
exponents, except in the ultimate vicinity of the critical point. In simulations, 
the chains in the blend are often too short to prevent critical fluctuations, 
and one finds Ising exponents. In the intermediate region between the weak 
segregation limit and the strong segregation limit, the contribution of loops 
competes with an increasing influence of chain ends. Interpolation schemes 



between the two limits have been sug geste d by Tang and Freed |101|| and more 
recently by Ermoshkin and Semenov[ |l02| . Numerically, the problem can be 
treated, e.g., within the self-consistent field theory [[2l"| p2| p3| . 



35 



From this discussion, it is clear that the properties of interfaces are strongly 
influenced by their microscopic structure. This interplay has attracted abiding 
theoretical pi ||, ||, |§ ||, [0n| and experimental |0§ |l0|, [Tl9| [l0|, [l06 
interest. 



07], 108, 109 1 interest. Apart from the factors already mentioned, packing 
and compressibility effects are again important, as well as chain orientations 
and chain-end distributions. These are all accessible to computer simulations. 
Unfortunately, simulation studies of polymer interfaces are computationally 
extremely demanding. Only recently have such simulations been carried out 
by different groups. The first is due to Reiter and coworkers and treats poly- 
disperse chains on a lattice model using a bond breaking algorithm) 11C ] . More 
recently, Ypma et al have performed simulations of a polymer interface con- 
fined between a hard wall and a free surface j70). Grest and coworkers study 
a polymer interface in a continuous model, in which demixing is driven by 
non-additive packing! Ill |. Similar to our simulations (8^, f75|, they also put 
particular emphasis on the analysis of the capillary- wave broadening. In gen- 
eral, the results from these simulations agree with our own results obtained 
with the bond-fluctuation model , which shall be discussed in some detail in 
the following sections. Simulations of copolymers embedded in a homopoly- 
mer interface are even more scarce. In many studies, the homopolymer inter- 
face has been approximated by an external potential[112, 113]. Simulations 



of a copolymer/homopolymer mixture have been performed by Pan et al[114|, 
however, the homopolymer concentration in their study is less than 10 %. 

We begin with the discussion of the interfacial properties of a homopolymer 
interface in a strictly symmetric binary blend. For a structurally symmetric 
polymer blend, the interfacial tension is accessible via the reweighting tech- 
nique. At the critical point, it vanishes with a critical power law of the form 
a = a(e — Ec)' 1 , where the critical exponent /i = 2v — 1.258 is given by the 
3D Ising universality class, and the critical amplitude takes the value a = 2.9 
for N = 32. In the strong segregation limit, the self-consistent field theory 
for Gaussian chains predicts in the long chain limit 



cssl = g^ 



( '2b \ + b A b B +b% 
bA + b B 



(55) 



where 6 a and bs denote the statistical segment length of the components. 
For finite chain length N, however, the interfacial tension is reduced by a 
factor 1 — a a -jy. There exist various analytical predictions for the constant 



= ln2|115|, a a = tt 2 /12[11§, a a = 1.35[|0l|, and a a = 21n2[|102| 



In Fig.|10| we present the ratio between the measured interfacial tension and 
the strong segregation estimate, and compare it to the self-consistent field 
theory [^|). The deviations from the strong segregation result are indeed well 
described by a correction of the proposed form, and the simulation data as 
well as the SCF results agree best with a constant 2 In 2. 

The influence of architectural asymmetry is illustrated in Fig.|ll|. Upon 
increasing the persistence length of the B component, the interfacial ten- 
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Figure 10: Temperature dependence of the interfacial tension for chain 
lengths N = 16, 32, 64. The ratio of the interfacial tension (measured us- 
ing the reweighting scheme) and the analytical prediction is presented to 
investigate chain end corrections. The lines correspond to self-consistent 
field calculations for compressible and incompressible systems. 
FromRef. ||. 
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sion increases [g3j. The arrow marks the result for symmetrical mixtures ob- 
tained via the reweighting technique. Note that packing effects, which are 
induced by the stiffness disparity, increase the incompatibility by a factor of 
Ax/x ~ 7 10~ 3 for / = Hence, they are not strong enough to explain 

the increase of the interfacial tension for this combination of chain length, 
cnthalpic repulsion and stiffness disparity. The simulation data in the fig- 
ure are compared to numerical self-consistent field calculations (SCF), which 
take due account of the chain conformations on all lengths scales, and to the 
analytical prediction for infinitely long Gaussian chains of different statisti- 
cal segment lengths (cf. Eqn. 55). The self-consistent field theory describes 



the dependence on the chain architecture almost quantitatively, whereas the 
analytical calculation captures the qualitative trend, but fails to predict the 
correct absolute values. 
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Figure 11: Stiffness dependence of the interfacial tension for chain length 
N = 32 for rather high incompatibility e = 0.05. Circles denote esti- 
mates via an analysis of the spectrum of interfacial fluctuations. The 
arrow marks the independently measured value (em via the reweight- 
ing scheme) for the symmetric blend. The data are compared to self- 
consistent field calculations which take due account of the chain archi- 
tecture. Also shown is the strong segregation prediction of the Gaus- 
sian chain model for blends with different statistical segment length by 
Helfand and Sapse (HS). The inset presents the relative effect of increas- 
ing the stiffness / of the B component. 
FromRef. ||. 
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These thermodynamic properties do not shed light on the microscopic 



structure of the interface[117 . However, factors like the width of the inter- 
facial region, the orientations of polymers on different length scales and the 
enrichment of vacancies or solvent influence the materials properties. Exper- 
iments indicate ]! 18| ], that entanglements in the interfacial zone are of major 
importance for the mechanical properties of the composite. The distribution 
of chain ends are important for the interdiffusion and healing properties at 
interfaces between long polymers. They also play an important role for reac- 
tions at interfaces. Chain ends are enriched in the center of interfaces, and 
this effects goes along with a depletion at a distance R g away from the mid- 
point of the profile. The shape of a polymer is a prolate ellipsoid, which is 
oriented in the presence of an interface. This is quantified by the orientational 
parameter 

where the outer index z at the brackets denotes the z coordinate of the vector's 
midpoint, and the inner indices the Cartesian components of the vector X. 
The apparent orientation profiles for the end-to-end vector R e and the bond 
vectors b are presented in Fig. Both the bond vectors and the end-to-end 
vector align parallel to the interface. Most notably, the orientation effects are 
the stronger the larger the length scale The Gaussian chain model does 
not predict any orientation of bond vectors, and the observed effect is rather 



small. The inset shows the predictions [52 of a worm-like chain model for the 



bond orientation. A rather small bending stiffness is sufficient to describe the 
simulation data (for chains with bending potential / — 0). On the end-to-end 
vector, however, the interface has a similar effect than a free surface or a hard 
wall, and induces pronounced alignment. 

In blends with stiffness disparity the persistence length of the semi- 
flexible component introduces a second microscopic length scale into the in- 
terfacial profile, which is independent from the interfacial width. In the weak 
segregation limit, the interfacial width is larger than the persistence length of 
the semi-flexible component and the interfacial behaviour can be described ap- 
propriately within the Gaussian chain model, when the stiffness dependence of 
the statistical segment length b(f) is accounted for. Upon increasing the stiff- 
ness disparity at fixed incompatibility, the interfacial width increases too. In 
the strong segregation limit, the interfacial widths becomes narrower than the 
persistence length and the local chain architecture becomes important. In this 
highly incompatible regime, the interfacial width decreases upon increasing 
the stiffness disparity at fixed x_/V, in qualitative contrast to the predictions 
of the Gaussian chain model. Self-consistent field calculations which take due 
account of the chain architecture on all length scales describe the simulation 
data almost quantitatively (5^]. Figure |f3| shows bond orientational profiles in 
a temperature region, where the stiffness disparity hardly affects the width of 
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Figure 12: (a): Orientation of bond vectors as a function of the dis- 
tance from the interface. Apparent profiles for chain length N = 32 are 
shown for various inverse temperatures e. The temperature ranges from 
the weak segregation limit to the strong segregation limit. The inset 
presents self-consistent field calculations in the framework of the worm- 
like chain model for different values of the bending rigidity r\ — 0.1,0.5 
and 1. and e = 0.1. (b): Apparent profile for the orientation of the 
chain's end-to-end vector as a function of the distance from the inter- 
face. Polymers align parallel to the interface. Note that the effect is 
much more pronounced than for the bond vectors. 
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the composition profile. Upon increasing the bending energy / of the semi- 
flexible component, the orientation of the semiflexible B polymers increases, 
whereas the flexible component remains unaffected. The width of the oricnta- 
tional profiles increases upon increasing /. The numerical self-consistent field 
calculations and the reduced interfacial profiles agree almost quantitatively. 

Moreover, we find enrichment of vacancies at the interface of our com- 
pressible blend. This segregation of vacancies can be qualitatively estimated 
as follows: By reducing the monomer density Sg, the system lowers the num- 
ber of unfavourable contacts between unlike monomer species. This reduces 
the excess energy density by ewz c Sg/2. The decrease of contact energy con- 
trasts the entropy loss due to the creation of a density fluctuation. The latter 
can be estimated as wSg 2 /2qkt, where kt denotes the isothermal compress- 
ibility. Balancing both contributions p3l| , we estimate the reduction of the 
total monomer density Sg to be of the order x^Tg/2. This simple argument 
is compatible with the simulation data for the total monomer concentration 
profile, and with self-consistent field calculations for compressible polymer 
blends. 

Finally in this section we discuss the profile of the number of intermolecu- 
lar contacts across the interface (Fig. |lj ]5l]] ). In the weak segregation limit, 
the number of contacts z c is independent of the distance from the interface, 
whereas it is strongly reduced in the center of the interface at high segrega- 
tion. This effect does not simply reflect the density reduction at the center 
of the interface, but is even more pronounced. It results from orientation and 
conformational changes of the macromolecules at the interface ]5l|: Polymers 
rearrange as to exchange energetically unfavourable intermolecular contacts 
with intramolecular contacts. The latter number is increased at the interface. 

4.3 Copolymers at interfaces 

We will now discuss the effect of adding copolymers to the interface. Such 
amphiphilic molecules are, of course, most likely to be found at the interface 
between A and B homopolymer regions, where the different blocks can ex- 
tend into the appropriate volume. This reduces not only their own enthalpy, 
but also that of the homopolymers which it displaces from the interface. The 
reduction of the interfacial tension has a strong influence on the morphol- 
ogy of the blend, whereas the entanglement of the individual blocks with the 
corresponding homopolymers increase the mechanical stability of the com- 
posite. The tendency to adsorb at the interface is balanced by the loss of 
entropy of mixing and of conformational entropy as the copolymers stretch 
to accommodate a greater areal density. Thus the adsorption at the interface 
competes with the creation of micelles in the bulk phases. Mean field theories 
predict that the addition of copolymers drives the system to compatibility in 
the weak segregation limit, while one encounters a first oder transition to a 
lamellar phase at strong segregation. In the intermediate regime, a complex 
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Figure 13: Orientation of bond vectors q&, for flexible A-polymers (solid 
lines) and semi-flexible _B-polymers (dashed lines) with bending energies 
/ = 0,0.5, 1, 1.5 and 2. The orientation of the i?-bonds and the length 
scale of ordering increases upon increasing /. (a) Monte Carlo results, 
(b) self-consistent field calculations The inset shows the self-consistent 
field results normalised to the maximum of q^. Note, that the width of 
the orientation profile grows upon increasing /. 
FromRef. @. 
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Figure 14: Apparent profiles of the intermolecular contacts for inverse 
temperatures e = 0.1, 0.07, 0.05, 0.035, 0.025, 0.019 and chain length N = 
32. Note the pronounced reduction of intermolecular contacts at the 
center of the interface in the strong segregation limit. 
FromRef. 
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phase behaviour is predicted by the theory. 

First, we discuss the adsorption of a small amount of copolymers at a 
strongly segregated homopolymer interface. The areal density v of copoly- 
mers is so small, that they do not overlap irvb 2 N/12 w 0.3. The apparent 
profiles of the individual segments are presented in Fig.[l5|. Chain ends of 
the homopolymers p e are enriched at the interface, whereas the density of 
middle segments is reduced compared to the total homopolymer density. The 
copolymers, however, exhibit the opposite trend. The middle segment, which 
joins the distinct blocks, is enriched at the center of the interface, whereas 
the chain ends stretch out into their corresponding bulk phases. The inset 
presents the results of the self-consistent field calculations, which qualitatively 
agree with the simulation data, except for the fact that they do not capture 
the capillary-wave broadening. The qualitative behaviour is also born out in 
experimental findings of Russell fll9| and co-workers. 

It is instructive to investigate the orientations of the copolymers at the 
interface. The normalised extensions of the end-to-end vector components 
parallel and perpendicular to the interface of the individual blocks are shown 
in Fig.|l6|. Right at the interface (z = 0) the individual blocks align parallel 
to the interface, similar to the behaviour of the homopolymers. In the minor- 
ity phase, the block shrinks, in order to reduce the number of unfavourable 
interchain contacts. Deep inside the majority phase, the block aligns perpen- 
dicular, because the other block pulls it towards the interface. On average, 
the orientation of individual blocks is however parallel to the interface. 

The orientation of the copolymer bonds is shown in Fig.[lj]. The link 
bonds between the different blocks are strongly oriented perpendicular to 
the interface. Their orientations increase as the distance from the center of 
the interface grows. The orientation of the bond vectors, however, rapidly 
decreases as we approach the ends of the blocks. The bonds in the middle of 
the blocks behave similar to the bond vectors in homopolymers. 

In sum, the shape of a copolymer resembles a dumbbell, the individual 
blocks are only mildly perturbed and resemble the homopolymer conforma- 
tions at the interface. The vector, which connects the centers of the distinct 
blocks, however, has a strong orientation perpendicular to the interface. 

The interplay between the interfacial thermodynamics and the phase be- 
haviour of the ternary homopolymer-copolymer mixture is best investigated 
in the semi-grandcanonical ensemble, which allows the simultaneous measure- 
ment of the interfacial free energy and excess quantities. Following Leibler[ 12C ] , 
we estimate the chemical potential of the copolymers in the bulk by: 

N8 i i = \np c + l - X N, (57) 

where the first term represents the translational entropy of the copolymer, 
and the second one the enthalpic repulsion between the homopolymer and 
the copolymer. The potential of a copolymer adsorbed at the interface, with 
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Figure 15: Copolymer segment density profiles for a symmetric diblock 
copolymer of length 32 at a highly segregated homopolymer interface 
N = 32 and e = 0.1. Apparent profiles are shown for the density of 
A and B monomers in the middle of the chain (pi/2, squares), at the 
end of the chain (p e , circles), and at one and three- fourths of the chain 
(pi/4, diamonds), and for the density of all copolymer monomers (p c , 
broken line). The inset shows the predictions of the self-consistent field 
theory for A monomers. Full lines show the predictions for segment 
density profiles P1/2, p e , an d P\/ 4, while the broken line presents the total 
density profile p c , and symbols compare with the MC results (notation 
like above). The distance from the interface is measured in units of the 
predicted width u>ssl in the strong segregation limit and densities are 
normalised appropriately. 
From Ref. fl7| . 
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Figure 16: Mean square end-to-end vector components of the copolymer 
blocks in their minority phase (A block in £?-rich phase or vice versa) 
and majority phase, in units of the average bulk value b 2 N/6, plotted as 
a function of the distance of the midpoint of the end-to-end vector from 
the interface. 
From Ref. |p!7 |. 



46 




z/w. 



SSL 



Figure 17: Orientation of bond vectors in a symmetric copolymer of 
length N — 32 at a highly segregated homopolymer interface (e = 0.1). 
The orientational order parameter q for end bonds, link bonds, which 
join the tw o blo cks, and neighbours of the link bonds are presented. 
From Ref. [pd7 |. 
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the A and B block extended into the appropriate homopolymer regions, is 
determined by the two-dimensional translational entropy of the joints at the 
interface, 

vN 

NS/ifuhx , (58) 

QWc 

where wc is the width over which the copolymer joints are confined to the 
interface, and v is the excess copolymer density at the interface. As shown 
in the segmental profiles (cf . Fig.[l5|) , the width is of the same order than the 



interfacial width. Semenov|121] has given the estimate wc = ttwssl/2 for 
strong segregation. In equilibrium, the chemical potential of the copolymers 
at the interface and in the bulk are equal, hence 

- ^ C exp(VY (59) 



N 1 \2 

The above considerations neglect the stretching penalty upon crowding copoly- 
mers at the interface. This becomes important when the copolymer concen- 
tration at the center of the interface is large, i.e., ^/ONv/gb ~ O(l). For 
higher interfacial excess, the number of copolymers does not increase linearly 
with the bulk concentration. The dependence of the interfacial enrichment 
on the bulk concentration for different temperatures is presented in Fig|lS| 



and compared with an analytical expression due to Scmcnov[121|. The in- 
terfacial excess increases upon adding copolymers in the bulk. At low bulk 
concentrations, we find rather good agreement with the theoretical descrip- 
tion, without any adjustable parameter. The saturation at higher copolymer 
concentrations is described only qualitatively. In this regime of intermediate 
segregation, the copolymers do not form a well defined stretched brush [fr5f, 
but rather an unstructured thick layer. Hence the theory overestimates the 
amount of stretching and the associated entropy loss. This goes along with 
an underestimation of the copolymer excess. 

The adsorption of copolymers also leads to a decrease of the interfacial 
tension a. This quantity is directly accessible via the probability distribution 
of the composition or the Gibbs adsorption isotherm 

vcNSfi, (60) 



knT 



where we have neglected the compressibility of the blend. The results of the di- 
rect measurement and the Gibbs adsorption isotherm are presented in Fig.[l9[ 
Both estimates agree nicely, indicating that the semi-grandcanonical simula- 
tion scheme establishes equilibrium between the bulk and the interface [f75f. 
From the slope at vanishing copolymer concentration, we can determine the 
effective width of the copolymer joint profile wc, which agrees rather well 
with Semenov's estimate wc = nw/2. In the weak segregation limit, the 
adsorption of copolymers leads to a compatibilisation of the blend. This 
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Figure 18: Segregation of symmetric copolymers of length TV = 32 at a 
homopolymer interface. The dependence of the excess areal density v of 
the copolymers on the bulk concentration is presented for various tem- 
peratures. The simulation data are compared without any adjustable 
parameter to Semenov's theory, which interpolates between the mush- 
room regime at low v and the brush-like regime for large excess densities. 
FromRef. |§. 
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second-order transition is indicated by the arrows in the figure. At higher 
segregation, we encounter a first order transition at which two homopolymer 
rich phases coexist with a spatially structured copolymer-rich phase. Both 
regimes are separated by a tricritical point J75|. The simulations indicate that 
the spatially structured phase is a swollen lamellar phase at high incompat- 
ibilities, and a microemulsion at intermediate incompatibilities. The latter 
has been established by inspection of the composition correlation functions, 
which are found to oscillate in space. 




Figure 19: Reduction of the interfacial tension upon increasing the 
copolymer's chemical potential Sfi. Symbols represent the value of the 
interfacial tension measured via the reweighting scheme, and lines give 
the result of the Gibbs isotherm, using the data from the previous fig- 
ure. The stability limits of the coexisting homopolymer-rich phases are 
indicated by arrows. Note that the addition of copolymers can reduce 
the interfacial tension by roughly two orders of magnitude, which is a 
typical value for amphiphilic systems forming microcmulsions. 
FromRef. |§. 

A detailed study of the bulk thermodynamics in the symmetric ternary 
homopolymer copolymer blend yields a tentative phase diagram over the 
whole range of copolymer concentrations and incompatibilities (see Figj2C|) 
It includes a two phase region at low copolymer concentration, where two 
homopolymer-rich phases coexist, and a disordered phase at low incompat- 
ibilities and high copolymer concentration. The disordered phase near the 
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tricritical point (TP) is a microemulsion. A swollen lamellar phase is found at 
higher incompatibilities and copolymer content. The first order phase transi- 
tion between the disordered phase and the lamellar phase is indicated by the 
double solid line. Unfortunately the transition between the microemulsion 
(D1S) and the (highly defective) lamellar morphology (LAM) could not be 
located accurately due to the finite simulation box. Hence, the double solid 
line is just meant schematic. However, the typical snapshots of the morphol- 
ogy also shown in the figure illustrate the qualitative difference between the 
lamellar phase and the region where a polymeric microemulsion is stable. 
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Figure 20: Tentative phase diagram of a symmetric homopolymer- 
copolymer mixture of chain length N = 32 on the isopleth (A/x = 0). 
The diagram includes the coexistence between two unstructured ho- 
mopolymer rich phases (2-phase) and a copolymer-rich disordered phase 
(D1S) at low incompatibilities. The disordered phase (DIS) is spatially 
structured (microemulsion) in the vicinity of the triple point (TP). At 
high incompatibility, we find a strong first order transition (with a large 
miscibility gap) between two homopolymer-rich phases and a lamellar 
copolymer-rich phase (LAM). The transition between the disordered 
phase (DIS, microemulsion) and the lamellar phase could not be located 
in the simulations, and the double solid lines are meant schematically. 
Typical snapshots of the A monomer density in a thin slice through the 
systems are also presented. For clarity the simulation box (L = 80) and 
three periodic images are shown. The arrows indicate the locations in 
the phase diagram. 



In the weak and strong segregation limit, the simulations agree with self- 
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consistent field calculations. At intermediate segregation, however, the nu- 
merical self-consistent field calculations by Janert and Schick Q predict the 
coexistence of two asymmetric highly swollen lamellar phases in the region, 
where the simulations find a polymeric microemulsion |75|] . Furthermore, the 
self-consistent field calculations predict a multi-critical Lifshitz point instead 
of a tricritical point. 

A discussion of the bending rigidity k of the interface provides a tentative 
explanation for the stability of the microemulsion observed in the simulations: 
At intermediate segregation, the interfacial tension a becomes very small in 
the vicinity of the triple line. Thus the bending rigidity controls the interfa- 
cial fluctuations, and t he p ersistence length £ p of the interface is of the order 
£ p ~ &exp(27rK/fcsT) [|l22| . The bending rigidity of the interface is accessible 
in the simulations via the spectrum of interfacial fluctuations |f75| (see Sec. 3). 
It is presented in Fig]2l]. The fluctuation spectrum for the pure homopoly- 
mer interface is well describable by the capillary- wave Hamiltonian ( |38| ) using 
the independently measured interfacial tension, and this description accounts 
also rather well for interfaces with a small amount of copolymer adsorbed. 
Upon increasing the copolymer excess, the interfacial tension decreases, and 
the deviations from the capillary-wave approximation ( |40| ) increase in turn. 
This indicates that the bending rigidity becomes important, and a satisfac- 
tory description of the simulation data is only possible with the full Helfrich 
Hamiltonian (p2[). In order to determine the bending rigidity k, we fit the 
simulation data to Eqn. ([l3|), fixing the interfacial tension to the value deter- 
mined via the reweighting scheme. Unfortunately, the bending rigidity cannot 
be extracted with a high accuracy, yet for our purpose here we only need an 
order of magnitude estimate. The values of the bending rigidity are small ( 
27tk/A;bT < 0.5) and increase with growing copolymer concentration. The 
order of magnitude and the way they depend on the copolymer co ncen tration 
agree qualitatively with analytical calculations on similar systems [ 123 1. 

The small value of the bending rigidity of a copolymer-laden interface sug- 
gests that interfacial fluctuations destroy the lamellar order predicted by the 
mean field theory, and lead to the formation of a microemulsion at intermedi- 
ate segregation. This int erpr etation is in accord with Monte Carlo simulations 
of Gompper and Kraus| 124 |, who found that the transition between a mi- 
croemulsion and a lamellar phase takes place at about 27tk/A;bT sa 16. Note 
that the chain length N = 32 corresponds to a degree of polymerisation of 
about 150 on an atomistic level. Our results are thus typical for rather short 
polymers, and the simulation data exhibit a behaviour between amphiphiles 
and high molecular weight polymers. 

On increasing the chain length N, the interfacial tension even decreases 
like 1/VjV at constant %iV, whereas the bending rigidity n is expected to 



grow[ 125 . Thus k dominates the fluctuation spectrum in the long chain length 
limit. De Gennes and Taupin ]l22] | argued that the lamellar phase becomes 
more favourable for long chain lengths. Therefore, the phase diagram should 
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Figure 21: Spectrum of interfacial fluctuations for a symmetric 
homopolymer-copolymer mixture at e = 0.054, i.e., intermediate seg- 
regation for three different chemical potentials of the copolymer (5/i = 
0.25, 1.375, 1.75 from top to bottom. The solid lines shows results of the 
capillary-wave approximation k — using the independently determined 
interfacial tension a. The arrows mark Semenov's estimate for the cut- 
off (?max — 2/w for Sfi — 1.375 and 1.75. The dashed lines correspond 
to fits according to the Hclfrich Hamiltonian, using the same interfacial 
tension as before. The fit values of the bending rigidity are indicated. 
FromRef. |§. 
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become more mean field like for longer chains. However, in the ultimate vicin- 
ity of the tricritical point, both the interfacial tension and the bending rigidity 
vanish and hence interfacial fluctuations are expected to become important 
for all chain lengths. Indeed, general arguments indicate that isotropic Lif- 
shitz points are unstable in three-dimensional systems ]126|. For very long 



chain lengths, Miiller and Schick) 75 have speculated that there might be a 
crossover from Lifshitz tricritical behaviour (£ p large, swollen lamellae) to or- 
dinary tricritical behaviour (£ p small, microemulsion) as one approaches the 
tricritical point. Recent experiments by Bates et al. (Tcfl found experimen- 
tal evidence for a bicontinuous polymeric microemulsion in an asymmetric 
homopolymer copolymer mixture in a region where the self-consistent field 
theory predicts an isotropic Lifshitz point. In fact, the experimental phase 
diagram of the copolymer homopolymer mixture closely resembles the phase 



diagram of small molecule mixtures of oil, water and amphiphilcs [ 127 



4.4 Thin films 

The behaviour of confined complex fluid is of practical importance for var- 
ious applications, e.g., adhesives, coatings, and lubrication. Confining sur- 
faces may alter the phase behaviour profoundly, and the interplay between 
interfaces separating coexisting bulk-like phases of a binary mixture and the 
confining walls has attracted longstanding theoretical interest JO], |Tij| . 

In the following, we discuss the behaviour of a binary polymer blend, which 
is confined between two hard, impenetrable walls. Neutral walls, which do not 
preferentially interact with any of the components, reduce the critical tem- 



perature of a symmetric polymer blend and likewise the miscibility gap[128| 



This has been investigated through extensive Monte Carlo simulations by 



Kumar |129[ and collaborators, and by Rouault[130| et al. The coexistence 
curve at the critical point is flattened compared to the bulk binodal, which 
indicates that the critical point belongs to the 2D Ising universality class. 
Upon approaching the critical point, the system crosses from mean field to 
Ising critical behaviour and passes from three-dimensional to two-dimensional 
critical behaviour. 

In general, one component may adsorb preferentially at the surface, such 
that the wall is coated with a layer of the component with the lower surface 
free energy. The structural and thermodynamical properties of these wetting 
layers are of practical importance and of fundamental interest in the statistical 
mechanics of condensed matter. At phase coexistence of the binary mixture, 
the surface free energy in the semi-infinite system undergoes a transition, at 



which the thickness of the adsorbed layer diverges. This wetting transition 15 
[l6| , |l7f may be continuous (second-order wetting) or the thickness may jump 
from a finite value to infinity at the wetting transition temperature. While 
the unmixing temperature grows linearly with the chain length, the wetting 
transition in symmetric blends is independent of the chain length. Unlike 
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the generic situation in mixtures of small molecules, wetting in polymeric 
systems thus occurs far below the critical point [ 131 1. The wetting behaviour 
in a binary polymer blend has been studied by Wang[132 and co-workers via 
Monte Carlo simulations in the semi-grandcanonical ensemble. 

In a thin film with non-neutral walls, two extreme cases are of special inter- 
est: In the asymmetric case, one surface attracts the A component, whereas 
the other wall prefers the B component. In the symmetric situation, both 
walls attract the same component with equal strength. We discuss their phase 
behaviour in turn: 

In the symmetric situation, the short range potential, which attracts pref- 
erentially the A component, shifts the coexistence chemical potential A/i away 
from its bulk value) 128]. At coexistence, the confined system phase separates 
laterally into A rich domains in coexistence with regions, in which A-rich lay- 
ers cover the surfaces, but the B-component prevails in the center of the film. 
In the temperature range between the critical temperature of the film and the 
wetting temperature, the thickness of the wetting layer at coexistence is deter- 
mined by the interplay of the repulsion between the wall and the AB- interface, 
which favours a thick wetting layer, and the shift of the coexistence chemi- 
cal potential, which suppresses the total amount of the A component in the 
film[133|. The semi-grandcanonical ensemble in junction with the reweighting 
methods permits an accurate location of the coexistence chemical potential, 
and a good characterisation of the phase behaviour over the whole tempera- 
ture range. The phase diagram in the bulk, and in a thin film is presented and 
compared to self-consistent field calculations in Fig. 22 83| . The simulations 
and the self-consistent field calculations agree qualitatively: The critical tem- 
perature is reduced, and at the critical point, the component favoured by the 
walls is enriched. We find a strong first order wetting transition far below the 
critical temperature T w w 0.2T C . Most notably, the A-pooi binodal is convex 
in an intermediate temperature regime. This curvature is the signature of the 
wetting transition in the semi-infinite system. 

In the as ymmetric case, the effect of the wetting transition is even more 
fundamental] 134, 135 1 : There always exists at least one interface which runs 
parallel to the walls and separates an A-rich phase from a S-rich one at 
the appropriate walls. In the temperature regime between the bulk critical 
temperature and the interface localisation-delocalisation temperature T C (D) 
of a film of width D, it is localised in the middle of the film. The total 
composition of this "soft-mode" phase fluctuates around 1/2. Hence, there is 
no symmetry breaking, and the system is in a one phase region. Below T C (D), 
two phases coexist in which the interface is localised close to one of the walls. 
The transition temperature T C (D) is smaller than the wetting temperature 
T w of the semi-infinite system, and approaches T w exponentially fast as the 
film thickness D grows. 

In the soft-mode phase, the deviations of the local interfacial position u 
from the middle of the film can be described by a capillary- wave Hamiltonian 
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Figure 22: Comparison between the Monte Carlo results (left figure) 
and the self-consistent field prediction (right diagram) for the phase be- 
haviour of a symmetric, binary polymer blend of chain length N = 32 
confined into a thin film of width D = 48 (MC) and D = 49 (SCF). 
Both surfaces attract the A component with equal strength e w — 0.16. 
The wetting temperature of the semi-infinite system is indicated by an 
arrow. The figures also include the bulk phase diagrams for comparison. 
From Ref. @. 



(cf. Eqn. |38|) augmented by the effective potential exerted by the walls onto 
the interface. We take the effective interaction between the interface and both 
walls, which binds the interface to the middle of the film, to be of the form 
V(u) = a(T - T c (D))exp{-\D/2)u 2 . Here 1/A denotes the effective range 
of the exponentially decaying potential between the interface and the wall, 
hence, exp(— AD/2) sets the strength of the interaction in the middle of the 
film. Therefore, the effective Hamiltonian reads 

H = J dx dy (|(Vm) 2 + a(T - T C (D)) exp(-AL>/2)M 2 | . (61) 

After Fourier decomposing the local interfacial position, we can calculate the 
spectrum of interfacial fluctuations 



L 2 (a(q) 2 ) k B T P U 



(62) 



As in the case of the free interfaces (Eqn. (41)), we integrate the spectrum of 



capillary waves from q m i n = 2it / L to g m ax, and obtain 



/ 2/ \\ k B T l + (g m ax£||) 2 

{u (x ' y)> = ^ ln Vi+(^„) 2 - (63) 

Note that the parallel correlation length £n plays a similar role than g m i n in 
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the previous analysis [^5|. It is given by 

fU = ^= = , eM+XD/A) (64) 



r^v y/2aa(T-T c (D)y 



Most notably, £m grows exponentially with the film thickness D. This is the 



characteristic signature of the soft mode phase) 134]. In a finite system. [p6[, £|| 
can exceed the lateral system dimensions. Depending on whether q m i n S,\\ = 
27r£||/L is much larger or much smaller than one, the spectrum is then cut off 
by the finite lateral system size L, or by the parallel correlation length £|| . 

The first limit, (? m m£|| ~^ °°j corresponds to interfacial fluctuations of a free 
interface. In the canonical ensemble, the constraint on the composition fixes 
the laterally averaged position u of the interface and Eqn. ( |63| ) can be used 
to determine the interfacial width. In the semi-grandcanonical ensemble, the 
average position of the interface is unconstrained, and one has an additional 
q = 0-mode, which corresponds to a translation of the interface as a whole, 
and which dominates the broadening. With the effective Hamiltonian (|Sl|), 
the fluctuations of the average interfacial position are then given by (u 2 ) — 
^bT^ 2 I (L 2 a). In the limit L -c £|| , the fluctuations are of the order D 2 itself, 
and thus w — D/2 in the semi-grandcanonical ensemble. 

In the second limit, <7mm£|| 0, the dependence of the parallel correlation 



length on the film thickness leads to an anomalous size dependence [137, 
of the apparent interfacial width. Using the convolution approximation ( |46] ) 
and the above equation for the lateral correlation length, we get 

2 2 ksT 2 ksTXD 
w « w + — — In g max £|i = w H — h const. (65) 

4(7 10(7 



Hence, the width of the apparent profiles increases like \/D for large lateral 
extensions. The dramatic dependence on the film thickness demonstrates the 
importance of careful assessment of interfacial fluctuations, when analysing 
apparent interfacial profiles. This finding is also relevant to interfaces in 
mixtures of of small molecules, and in particular to the interpretation of ex- 



perimental data[138|. In the symmetric case (capillary condensation) or in 
the presence of long ranged (van der Waals) interactions, the dependence of 
the interfacial width on the film width D is, however, only logarithmic. 

The structure and profiles across the film in the soft-mode phase have 
been studied with extensive Monte Carlo simulations by Werner et al.fjof. 
The system considered is a symmetric, binary blend of polymers with chain 
length Na = Nb = 32, confined to a thin film of width D. The left surface at- 
tracts A monomers, which are less than 2 lattice spacings away, and repels B 
monomers, with an interaction of strength e w = 0.1. Likewise, the right sur- 
face prefers B monomers and repels A segments. The temperature 1/e = 33.3 
was chosen higher than the wetting temperature in the semi-infinite geometry, 
hence the system is without doubt in the soft-mode phase. In Fig. E3l the 
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apparent intcrfacial width w is shown as a function of the film thickness D for 
very large lateral system size L — 256 in the semi-grandcanonical ensemble. 
The large lateral system extension ensures, that the fluctuation spectrum is 
cut off by the parallel correlation length £|| . The Monte Carlo data are com- 
pared to the prediction (|65|), where the r ange of t he interfacial interaction A 



is taken to be 1/A = £(1 + fcsT/87r^ 2 cr)||l35j, 136 1 with the bulk correlation 
length £, and the constant term is taken to be zero. Both, the bulk corre- 
lation length and the interfacial tension have been measured independently. 
The intrinsic width has been extracted from a self-consistent field calculation 
in an infinite system j52|. The Monte Carlo data and the predictions agree 
nicely for thick films; however, there are deviations for very thin films. This 
can be understood from the following consideration: If the film width D is not 
very much larger than the radius of gyration R g (N — 32) 7, the intrinsic 
profile po is squeezed. Such an effect has been observed in simulations [|S6| 
and in self-consistent field calculations ||| . It can be analysed via a block 
analysis (see Sec. 3) by inspection of the width of the reduced profiles as a 
function of the block size B. If B is larger than the short wavelength cut-off 
of the interfacial fluctuation spectrum l/(7 ma x, the reduced profiles are broad- 
ened by capillary waves. One can attempt to identify an effective intrinsic 
width with the value wo(D) attained at the subsystem size B = 8 « R g . As 
shown in Fig. |2^, Eqn.(p5|) in combination with this effective width wq(D) 
quantitatively accounts for the deviations at small film thicknesses. 

The relationship between the width of the apparent interfacial profiles 
and the film thickness follows from general considerations on the statistical 
mechanics of interfaces, and is not restricted to polymeric systems. Similar ef- 
fects have been observed previously in simulations of simple Ising models[137|. 
However, polymer blends are particularly suited for investigating phenomeno- 
logical concepts of surface enrichment, wetting, interfacial localisation-delocalisation 
transitions and capillary condens atio n, because the chain length N constitutes 



an additional control parameter [131 . On increasing the chain length N, one 
can reduce bulk composition fluctuations, which are ignored in most phe- 
nomenological approaches. There exist powerful analytical tools to describe 
the bulk and interfacial behaviour in the long chain length limit. In addition, 
the phenomena occur on larger length scales due to the large size of the poly- 
mer coils, which facilitates applications of several experimental techniques. 
In fact, the anomalous size dependence of the apparent interfacial width has 



first been observed experimentally in a polymer system by Kerle et al.[138 



5 Conclusions and outlook 

We have reviewed extensive computer simulations on the bulk thermodynam- 
ics and interfacial structure of polymer blends in the melt state. Simulations 
of polymer blends and interfaces are considerably more exacting in computa- 
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Figure 23: Monte Carlo results for the apparent interfacial profile width 
of a symmetric, binary polymer blend N = 32 in a thin film of thickness 
D. The confining surfaces attract different components of the blend 
with an interaction potential e w = 0.1. The temperature e = 0.03 is 
well above the wetting temperature. Shown is the squared interfacial 
width w 2 as a function of the film thickness D for large lateral system 
extension L = 256. The straight line shows the expected a capillary- 
wave broadening, when the interfacial fluctuation spectrum is cut off by 
the thickness dependent parallel correlation length £y . 
FromRef. @. 
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Figure 24: Same data as in the previous figure, but the squeezing of the 
intrinsic interfacial width due to the confining surfaces is accounted for 
by subtracting the estimated intrinsic width Wq(D). The straight line 
corresponds to the expected capillary-wave broadening. 
FromRef. |§. 
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tional terms than those of small molecular fluids or magnetic systems. The 
difficulties stem from the difficulty of dealing with the widely spread time 
and length scales caused by the extended structure of the macromolecules. 
The accurate determination of the macroscopic behaviour while retaining the 
detailed atomistic chemical structure is not feasible even with state-of-the- 
art supercomputers. Yet, by a careful choice of the statistical ensemble, use 
of recently developed simulation techniques and sophisticated data analysis 
methods, simulations of coarse-grained polymer models yield valuable insight 
into the structural and thermodynamic properties of polymeric composites. 

Simulations contribute to the identification of important effective param- 
eters, which determine qualitatively the bulk and interfacial properties. They 
can examine the influence of the coarse-grained chain architecture on the 
phase and interfacial behaviour. As an example, we have discussed how stiff- 
ness disparity |Q between the components gives rise to a positive contribution 
to the Flory-Huggins parameter \- This finding is in agreement with field 
theoretical considerations pj[ and P-R1SM calculations |37j . Furthermore, we 
have shown that the persistence length affects the interfacial behaviour in 
highly incompatible blends. For weak incompatibilities, the effect of the stiff- 
ness can be modelled within the Gaussian chain model. This model incor- 
porates two parameters, the incompatibility %iV, and the length scale R set 
by the polymer radius of gyration. In the strong segregation limit, the width 
of the interface scales as R/\/xN. However, for large incompatibilities the 
interfacial width becomes of the order of the persistence length. On this 
length scale, the conformations are not Gaussian, but the bond vectors are 
strongly correlated. In this regime simulations and self-consistent field calcu- 
lations show that the interfacial width decreases upon increasing the stiffness 
disparity, in qualitative contrast p3| to the predictions of the Gaussian chain 
model. 

Many static and dynamical properties on various length scales are simulta- 
neously accessible in simulations. Therefore the simulations provide a detailed 
picture of the interfacial properties in binary and ternary polymer mixtures. 
Both structural informations (e.g., the orientation of bond vectors, or of whole 



chains [fell p2, 117], the adsorbtion of copolymers |75|, the enrichment of vacan- 



cies) as well as thermodynamic properties (e.g., the interfacial tension or the 
bending rigidity) are obtained simultaneously and hence permit a quantitative 
comparison with other theoretical approaches like scaling descriptions 111, self- 



consistent field calculatio ns |21| , [22| , |23|, |24j] , density functional theories |25|, [101 
and P-RISM calculations 



Simulations of the interfacial behaviour of symmetric block copolymers [1 17 1 
yield a detailed picture of the structure. If the areal density of copoly- 
mers adsorbed at the homopolymer interface is small, the copolymers assume 
dumbbell-like conformations. The individual blocks extend into the appro- 
priate bulk phases and resemble only mildly perturbed homopolymers, in 
particular each block is aligned parallel to the interface. This finding is in 
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qualitative agreement with available experimental data by Russell et al.[119|. 

In general, simulations and mean field treatments agree qualitatively, be- 
cause the extended structure of the polymers suppresses fluctuations, and 
microscopic details can be qualitatively described with a few effective, coarse- 
grained parameters. The simulation can access the way in which the mean 
field behaviour is approached upon increasing the chain length |j49f . However, 
one important source of discrepancy between simulations and mean field the- 
ories are interfacial fluctuations and we have highlightened two consequences: 

In ternary homopolymer copolymer blends |75[| , interfacial fluctuations de- 
stroy the order in highly swollen lamellar phases and lead to the fluctuation- 
induced formation of a microemulsion. Simulations find tentative evidence 
for a tricritical point, where the mean field treatment [jl2[ predicts a Lifshitz 
tricritical point in symmetric systems. Upon increasing the chain length, 
however, there may be a crossover from Lifshitz criticality to ordinary tri- 
critical behavioural. This "catastrophic failure of mean-field theory" Jl0| is 
also observed in recent experimental studies by Bates and coworkers. These 
experiments on an asymmetric ternary blend find a bicontinuous microemul- 
sion close to the region, where the mean-field treatment predicts an isotropic 
Lifshitz point. The mesoscopic structured yet macroscopic homogeneous com- 
posites possibly have unique and useful mechanical or electrical materials 
properties Q. 

Interfacial fluctuations profoundly influence apparent profiles measured in 
simulations [^6). They lead to a pronounced broadening of the width of the 
measured profiles. In a thin film with asymmetric surfaces attracting different 



components, there exists a temperature regime[134 , in which one interface 



separating an A rich domain and a B rich one runs parallel to the surfaces 
in the middle of the film. Effective interactions between the interface and 
both walls pin this interface to the middle of the film, and generate a parallel 
correlation length £11 which grows exponentially with the film thickness D 
in in the case of short range interactions. This large length scale acts as 
a cut-off for the spectrum of capillary waves, and consequently limits the 
broadening. The apparent interfacial width increases like \JD. If we were 
to include long range van der Waals interaction in our simulations, we would 
expect a logarithmic dependence of the apparent interfacial width on the 
film thickness. The effe c t is also pertinent to the interpretation of recent 
experimental studies[138, 139, 14C ] . 

In view of these findings, it would be of interest in future simulation stud- 
ies to investigate the role of the short length scale cut-off of the interfacial 
fluctuations spectrum. A study of the way how the capillary- wave broadening 
depends on the chain length, and on the incompatibility, might help to solve 
the question, whether and how the length scales of the lateral fluctuations 
contributing to the intrinsic profile or to the capillary-wave broadening can 
be separated from each other. This might contribute in general to the in- 
terpretation of intrinsic profiles calculated in analytical approaches. Another, 
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more subtle point, concerns the way the shape of the profiles is affected by the 
capillary- wave broadening. The intrinsic profiles obtained, e.g., with a block 
analysis, should be reasonably well described by a tanh profile, at least in the 
weak segregation limit. Preliminary results indicate that the intrinsic profiles 
are indeed closer to a ta nh s hape than to an error function shape as the crit- 
ical point is approached) 14C]. The shape of the apparent profiles is given by 
the convolution of the intrinsic profile and the gaussian curve associated with 
the capillary broadening, which is closer to an error function in large enough 
systems. The latter error function shape has been established for t he ca se of 
capillary-wave broadening in Ising models only recently by Moseley[141 . 

Moreover, the dynamical evolution of interfacial profiles, after two flat 
interfaces have been brought into contact, is relevant in some experimental 
situations. Simulations might investigate whether there is a separation of time 
scales on which an intrinsic profile develops which is then, on a different time 
scale, broadened by capillary waves. 

A few other interesting questions shall be sketched briefly: The interfacial 
properties of random or branched copolymers have attracted considerable in- 
terest because of the extensive commercial applications. Recent interest stems 
from experimental progress of the reactive polymer blending techniques[i42|. 
This has been subsequently studied theoretically! 143, 144, 145, [146| for rather 
idealised situations. Other important problems are related to the single chain 
dynamics in blends. In the miscible case, e.g., simulations could test recent 
theoretical predictions by Schweizer[147| and co-workers on the coupling be- 
tween the thermodynamics and the single chain dynamics. 

Furthermore, recent attempts to construct even more coarse-grained model |l48| 
which incorporate the instantaneous shape of a polymer only via a continuous 
distribution of its monomer density, might offer promising routes to investi- 
gate the large length scale properties of macromolecular composites. This 
might shed some light onto the mechanical properties of polymeric interfaces, 
which are of long standing practical and experimental interest. 
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